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This  document  describes  the  Airport  Vicinity  Air  Pollution  (AVAF) 
Model  modified  by  the  Energy  and  Environmental  Systems  Division  of  Argonne 
National  Laboratory  for  the  Federal  Aviation  Administration  under  Inter- 
Agency  Agreement  DOT-FA71WA1-223 . This  version  of  the  model,  Airport 
Vicinity  Air  Pollution  hod  el  Abbreviated  Version,  provides  a 'tflret- 
guess*'  estimate  of  an  airport  air  pollution  distribution  based  on  more 
general  Input  Information  than  Is  described  In  Report  No.  FAA-RD-75-230, 
the  model's  primary  form.  ' 

The  user  must  establish  a right-handed  Cartesian  coordinate  system 
as  a basis  for  describing  an  airport  layout,  then  decide  whether  to  use 
supplied  default  values  applicable  to  the  types  of  altcraft  and  airport 
fac lilt les . 

' The  estimate  of  pollutant  concentrations  are  computed  for  selected 
Individual  receptor  locations  based  on  the  dispersion  of  pollutant  emissions 
from  area  and  line  type  sources.  Pollutant  emission  levels  due  to  each 
member  of  a source  type  are  tabulated.  The  results  of  a one  one-hour 
modeling  period  of  aircraft  arrival-departure  activities,  airport  ground- 
support  services,  and  airport  access  vehicle  roadways  are  tabulated  for 
each  receptor  location  and  pollutant  apecles. 

''The  document  presents  a discussion  of  the  theoretical  considerations 
fundamental  to  the  AVAP  Model  Abbreviated  Version  as  well  as  program 
flow  diagrams  essential  to  understanding  the  theory.  Required  sequential 
card  data  Input  to  the  program  Is  Illustrated  along  with  substitution 
options  for  program  constants.  Additionally,  an  example  problem  and  a 
program  listing  are  provided. 
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1 INTRODUCTION 


The  Airport  Air  Pollution  Model,  Abbreviated  Version  was  developed 
at  Argonne  National  Laboratory  and  executed  on  the  IBM  360  Model  195  within  a 
core  region  of  250K  bytes  while  utilizing  less  than  one  minute  of  computer 
time.  The  FORTRAN  coded  program  is  yet  another  version  of  the  Airport 
Air  Pollution  (AVAP)  Model^  modified  to  provide  a first-guess  estimate 
of  an  airport  air  pollution  distribution. 

This  document  discusses  the  theoretical  considerations  fundamental 
to  the  model,  describes  data  incut  reauirements  for  usine  the  Droeram,  and 
illustrates  the  use  of  model  by  presenting  an  example  problem.  Therefore, 
users  of  the  document  should  be  technical  personnel  concerned  with  using  the 
computer  code  for  assessing  air  pollution  impacts  of  commercial  and  general 
aviation  airports.  For  more  comprehensive  studies,  the  model's  primary  form 
should  be  used. 

In  Section  2,  the  collection  of  formulas  that  form  the  basis  of  the 
computer  program  are  presented.  Detailed  descriptions  are  to  be  found  for 
area  source  and  line  source  dispersion  formulas.  Although  they  are  not 
used  in  this  abbreviated  version  of  AVAP,  the  Gaussian  plume  formula  for  a 
point  source,  the  Carson-Moses  plume  rise  formula,  and  the  Briggs  formula 
for  estimating  stack  downwash  effects  are  presented. 

Section  3 provides  guidelines  for  preparation  of  the  input  data. 

The  formats  and  meanings  of  the  input  parameters  are  given  along  with  their 
relative  card  position  in  the  input  sequence.  Almost  all  of  the  input 
parameters  that  appear  in  this  section  havedefault  values.  Those  that 

do  not  have  default  values  are  marked  in  Table  3.1.  with  the  asterisk  (*) 
character. 

Section  4 provides  a description  and  a table  of  the  default  values 
of  the  input  parameters.  All  of  the  default  values  appear  in  the  BLOCK  DATA 
subprogram  of  the  computer  program;  all  or  some  of  them  can  be  changed  by 
the  appropriate  data-card  Insertion. 
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Section  5 presents  a complete  example  problem  that  details  the 
resultant  output  from  the  model.  Appendix  A illustrates  the  conversion  of 
pollutant  concentration  units  and  Appendix  B gives  the  program  flow 
diagrams  of  selected  computational  modules. 

A complete  listing  of  the  FORTRAN  coded  program  is  found  in 
Appendix  C. 
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2 FORMULAS  FOR  CONCENTRATION  CALCULATIONS 

The  formulas  presented  in  this  section  will  apply  to  both  forms  of  the 
AVAR  Model  except  that  Point  Source  calculations  are  not  required  when  using 
the  abbreviated  version.  They  form  an  outline  of  the  most  basic  set  of  equations 
for  concentration  calculations  but  are  not  to  be  regarded  as  the  complete  set 
that  covers  all  practical  conditions  of  application.  The  user  of  the  AVAP  Model 
should  consult  the  references  quoted  in  this  document  for  more  complete  descrip- 
t i ons . 

2.1  DISPERSION  COEFFICIENTS  0 and  0 

y z 

The  dispersion  coefficients  a and  O (miles)  indicate  the  amount  the 

y z 

pollutant  plume  has  spread  (dispersed)  after  leaving  its  source.  To  avoid  un- 
realistic behavior  of  the  o functions  at  very  high  and  very  low  wind  speeds, 
the  following  formulas  are  used: 

°y  = oy(T)  = Max  OyT(T),Oyx(x)]I  (1) 

and  a = 0 (T)  - Max  !a  _(T)  ,0  (x)  , (2) 

Z 2 Z 1 ZX 

where  x is  the  downwind  distance  and  T = x/u  is  the  travel  time.  In  Eqs.  (1) 

and  (2),  a _(T)  and  a (T)  denote  the  travel-time-dependent  dispersion  coeffi- 
yT  zT 

cients.  Curves  of  ayT(T)  for  different  Turner  stability  classes  are  displayed 

in  Figure  2.1.  Turner's  original  values  for  two-hour  sampling  time  have  been 

0 . 2 

converted  to  one-hour  sampling  time  by  multiplication  with  the  factor  (1/2) 

= 0.87.  No  conversion  factor  for  sampling  time  is  applied  to  the  OzT(T)  values 
which  are  plotted  in  Figure  2.2.  Figures  2.3  and  2.4  show  the  downwind-distance- 
dependent  dispersion  coefficients  0yx(x)  and  0zx(x)  derived  from  curves  in  Tur- 
ner's Workbook  for  10-minute  sampling  time,  by  multiplying  the  original  o (x) 

2 2 
values  by  (60/10)'  and  the  original  02X(X)  values  by  (20/10)’  . In  doing  so 

we  have  assumed  that  the  vertical  dispersion  coefficient  is  insensitive  to 

sampling  times  beyond  20  minutes,  as  suggested  by  Slade. 

2.2  POINT  SOURCE  DISPERSION  EQUATIONS 

The  short-term  average  concentration  at  the  recptor  point  (x,v,z)  due 
to  a point  source  at 


AxS  U‘i  • duO 


MILES 


Figure  2.4.  Itownwind-Distance-Dependent.  Vertical  Dispersion 
Coefficients  for  1-Hour  Sampling  Time. 
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y * crosswind  distance  from  plume  centerline, 

Z - distance  above  ground, 

u - wind  speed  at  source  height  H (See  Sec.  2.5). 
q = emission  rate  of  the  point  source  in  mass  per  unit  time, 
o = horizontal  dispersion  coefficient, 

y 

a =>  vertical  dispersion  coefficient 
z 

In  the  AVAP  Model,  plumes  emitted  from  point  sources  are  assumed  to  have 


finite  initial  dimensions  (0yQ  - 0 
downwash  does  not  appear). 


100  ft.  at  the  source  point  (0,0,H)  when 


If  there  is  no  interference  from  the  ski  lid,  (see  Sec.  2.6  for 
its  effects) the  source  height  is  generally  computed  as  the  sum  of  the  physical 
stack  height  h and  plume  rise  Ah  given  by  the  Carson-Moses  formula:5 


Ah  = A(Qh)1/2/Us(ft), 


s » wind  speed  (mi/hr)  at  physical  stack  height, 

A » 0.870  for  stability  classes  1 to  3, 

sgl  - 0.354  for  class  4, 
sp  - 0.222  for  class  5,  and 

Qh  - heat  emission  rate  of  the  point  source  in  Btu/hr. 

Stack  downwash  is  assumed  to  occur  if  Us  ’ 6 mph.  The  physical  stack 
height  h is  then  replaced  by  Briggs'  modified  stack  height:** 


Zs  + 2 <VU.  " 1-5)Ds* 

Zg  lfVs>1.5Us. 


h'  “ modified  stack  height  after  stack  downwash, 

- stack  effluent  velocity,  and 
D»  “ stack  diameter. 

Alsu,  with  stack  downwash  the  initial  plume  dimensions  are  reassigned  in  the  AVAP 
Model  io  simulate  average  city  block  size  and  building  height  (oyo  - 250  ft; 

*>,,  - 40  ft). 
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For  point  sources  with  poor  aerodynamic  characteristics,  such  as  vents 

or  very  short  stacks  on  buildings,  o Is  automatically  assumed  to  be  2‘>0  ft 

and  the  building  height  h Is  used  to  compute  a (a  ■ h/1.2).  If  h Is  not 

zo  zo 

supplied  by  the  user  the  default  value  of  40  ft  Is  used  for  <3ro. 

2.1  AREA  SOURCE  DISPERSION  EQUATIONS 

In  the  AVAP  Model,  two  sets  of  dispersion  equations  are  used  for  area 
sources.  For  convenience,  urea  sources  are  classified  as  "near"  sources  and 
"far"  sources  depending  on  the  relative  location  of  the  receptor  and  the  area 
source . 

First,  the  critical  distance  for  mixing,  x^,  is  computed  from  the  mixing 
height  value  L by  the  equation 

0 (T  ■ x /u)  ■ 0.47  L.  (6) 

z c 

If  the  critical  distance  measured  upwind  from  the  receptor  Is  downwind  of  the 
downwind  edge  of  the  area  source,  the  area  source  Is  defined  as  a "far" 
source.  If  the  critical  distance  measured  upwind  from  the  receptor  extends 
beyond  the  upwind  edge  of  the  source,  the  entire  area  source  Is  treated 
as  "near."  There  will  be  cases  In  which  an  area  source  Is  partitioned 
Into  a "near"  and  a "far"  source  relative  to  the  receptor.  The  detailed  logic 
for  area  source  classification  ts  presented  In  Appendix  B on  p.  B-14. 

Tlie  far  area  sources  are  treated  conventionally,  with  the  horizontal 

and  vertical  dispersions  represented  separately,  each  by  an  upwind  virtual  point 

source,  so  that  at  the  downwind  edge  of  the  area  source  and  at  the 

center  cl  « . For  an  area  source  of  side  length  d and  vertical  spread  h, 

z zo 

the  initial  Gaussian  widths  are  given  by  Oyu  - d/2.4  and  Ozo  - h/1.2. 

For  a complete  discussion  of  the  treatment  of  near  area  sources,  the 
reader  Is  referred  to  a separate  report.^  Briefly,  it  Is  to  treat  the  area 

source  In  such  a manner  that  the  z component  Is  represented  by  an  upwind  virtual 
line  segment  along  the  wind  direction  Instead  of  u single  point,  with  the 
simplifying  assumption  that  Oy  is  held  constant  over  the  area  source.  The 
/ - component  is  analytically  integrated  over  the  line  segment.  Therefore, 
for  receptor  locations  Immediately  downwind  of  the  area  source,  the  model  Is 
expected  to  give  more  realistic  concentration  profiles  than  the  conventional 
area  source  model. 
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The  short-term  average  concentration  at  tl»«»  receptor  point  (x,y,t), 

a 

due  to  an  area  source  having  its  geometrical  center  at  (0,0,2*),  la  given  by 


vy 

Xa  “ 2it  cr  d(l-b)  cxp 

y 


exp 


1 


where , 


F - 


<VT2> 


1 fl 

2 I 0 

V y 

“,<V  ! 


W 


+ exp 


, i Z + Z 

1 i a 

2 \ITJt2) 


(7) 


en  (0  (T  )/cr  (T  )) 

b VVV 

If  b - 1,  the  factor  F/(l-b)  in  Kq.  7 should  be  replaced  by  Tj  ^(Tj/Tj)/ 

a (T,). 
z 1 

In  the  above  equations,  denotes  the  sum  of  (pseudo  travel  time 
corresponding  to  d ) and  the  travel  time  from  the  downwind  edge  of  the  area 
source  to  the  receptor,  and  T.,  denotes  the  sum  of  and  the  travel  time  from 
the  upwind  edge  of  the  area  source  to  the  receptor.  If  the  receptor  Is  inside  the 
area  source,  T « 0 . is  computed  with  a travel  tlmoqo  that  Oy  • d/2.4  at. 

tho  Jwnuind  od»«  of  tho  sourco.  or  It  th.  receptor  U lnoldP  th.  soorco  It  t» 
assumed  to  be  d/2.4  at  the  receptor. 


2.4  FINITE  LINK  SOURCE  DISPERSION  EQUATIONS 

The  short-term  average  concentration  at  the  receptor  point  (x,y,z) 
due  to  a finite  line  source  with  inclination  angle  0 relative  to  ground, 
axlmuthal  angle  ^ relative  to  the  wind  vector  anil  its  end  points  at  (0,0,0) 
and  (l.cos  0 cos  <b , Lcos  0 sin  i|>,  Lsin  (0  is  given  by 

Xo  " _ 


.J, 


(8) 


2tt  uo  o 
y z 


with 


2 I 

- (A/2)  .A.^  n|oxp  (b^-C^  ' erf^Bt+L/A)  - orf  ( B()| 


2 2-1  2-t 

A • /2  it  a cos  0 sin*'  J>a  “ + sin^O  5 z 
v z z y 


A-  1/2 


J 


wliero 


1 \ 


r 


cost)  slnij>  O ‘ + (Z-H)  slot)  a 2 
Z y 

> 

- 2 - 2 
cost)  slnt  0 + (Z+H)  sinO  a 

z V 


The  above  formula  is  applicable  only  under  certain  conditions.  For  a 
complete  discussion  of  the  various  criteria  and  the  reasons  behind  them,  the 
user  is  referred  to  Kefs.  8 and  9,  Normally,  when  the  line  source  subtends  a 
sufficiently  large  angle  relative  to  the  wind  vector,  for  example,  a uniform, 
horizontal  line  with  4>  greater  than  45°,  the  formula  Is  used  without  the  seg- 
mentation scheme  discussed  in  the  references  mentioned  above.  The  dispersion 
coefficients  0 and  o are  evaluated  with  an  effective  downwind  distance  cor- 

y z 

responding  to  a point  on  the  line  that  is  directly  upwind  of  the  receptor. 
When  the  relative  angle  is  very  small  (<  10°),  the  following  approximation 
formula  is  used: 


Again,  for  a long  line  source,  the  line  is  divided  into  shorter  segments,  and  an 

effective  J and  a for  each  segment  la  evaluated  for  the  downwind  distance 

y z o o 

corresponding  to  the  midpoint  of  the  segment.  For  angles  between  10  and  41 

Eq.  8 is  used  with  the  sgmentation  scheme. 

2.1  WIND  PROFILE  LAW 

To  convert  the  wind  speed  measured  by  anemometer  at  a local  airport 
(typical  Iv  30  ft  above  ground)  to  that  at  the  physical  stack  height  h or  to 
that  at  the  effective  source  height  H we  use  a power  law  relation  of  the  form: 

u ( Z ) - u(Z  ) (?./?.  )P.  (»°> 

o o 


/ _ - \ -2 
B - - 2 O O 1 
1 ' y z 


B2  “ "2  ( 0 


_ •> 


° ) 


V z 


+ iZzil 1 
20  2 


+ I2±H 1 


.A. 

.A. 

2 1/2 

I 

1/2 


i y 


20 


- 2 
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The  exponent  P,  as  determined  by  DeMarrais  , depends  on  the  stability  class 
and  is  given  In  Table  2,1. 

Table  2.1  Exponents  for  Wind  Profile 


Stability  Class  P 


Unstab le 

(<3) 

0.2 

Neutral 

(4) 

0.3 

Stable 

(5) 

0.4 

Very  Stable 

(6) 

0.5 

J . 6 EFFECTS  OF  SKY  LID 


The  lower  surface  of  an  elevated  inversion  layer  is  referred  to  as  the 
sky  lid,  and  as  a rough  approximation  it  is  assumed  to  act  as  a perfect  reflec- 
tor. Thus  the  plume  would  generally  be  reflected  repeatedly  from  two  parallel 
surfaces,  the  ground  and  the  sky  lid.  If  the  reflected  plumes  from  a point 
source  are  represented  by  multiple  image  sources  (see  Fig.  2.5),  then  the  net 
concentration  at  the  receptor  can  be  expressed  in  terms  of  a multiple  image 
series : 


N 

■ Z 

1-0 


V 


(ID 


where  \ - concentration  due  to  original  source,  and 

P° 

- concentration  due  to  I1*'  image  source. 

The  above  formula  Is  used  with  N - 6 (i.e.,  three  Image  sources  below  ground 

and  three  above  the  sky  lid)  when  the  downwind  distance  x of  the  recoptor  from 

the  source  is  less  than  2 x as  calculated  from  Eq.(6).  For  downwind  distances 

x ' 2 x , the  net  concentration  is  calculated  with  the  assumption  of  full  mixing 
- c 

in  the  mixing  layer,  and 


3 

/2n  Lu  0 


(12) 


..  i 
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For  some  point  sources  (e.g.,  power  plants  with  tall  stacks),  the 
calculated  effective  emission  height  H may  be  greater  than  the  lid  height, 
especially  when  the  latter  is  small.  These  cases  are  not  eliminated  from  con- 
sideration in  the  AVAF  Model,  but  are  simulated  by  using  dispersion  coefficients 
and  plume  rise  for  stably  stratified  air  (stability  class  5).  Depending  on 
the  relative  position  of  the  lid,  the  physical  stack  height, and  the  computed 
effective  stack  height,  the  model  will  in  some  cases  reassign  the  effective 
stack  height  or  the  lid  height  values. 

For  near  area  sources,  there  are  no  sky  lid  effects  by  definition. 

For  far  area  sources,  the  sky  lid  effects  are  handled  the  same  way  as  the 

point  sources. 

The  sky  lid  effects  on  line  source  dispersion  are  handled  by  a simple 

linear  interpolation  scheme.  For  effective  downwind  distances  (see  Section  2.4 

for  definition)  less  than  the  critical  distance  xfi  , no  sky  lid  effects  are 

assumed;  and  for  effective  downwind  distances  greater  than  2 xc» uniform  vertical 

mixing  is  assumed  by  using  the  integrated  expression  for  the  Z-component* 

For  downwind  distances  between  x ^ and  2 x£  a linear  interpolation  is  applied 

to  the  concentration  values  obtained  at  x and  2 x . 

c c 
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3 DESCRIPTION  OF  INPUT 

The  first  six  (6)  cards  in  the  data  deck  (see  Table  3.1)  are  the 
problem  title  and  program  parameter  input  cards  required  for  every  appli- 
cation of  the  model. 

Card  1 is  a single  card  of  input  having  no  computational  value  in- 
ternal to  the  code.  It  may  be  punched  or  left  blank.  If  punched,  the  same 
will  appear  as  the  first  line  of  printed  output. 

Card  2 provides  at  least  one  identifying  symbol  chosen  from  the 
following  set: 

CO  - Carbon  Monoxide 

THC  - Total  Hydrocarbons 

NOX  - Nitrogen  Oxides 

PART  - Total  Suspended  Particulate 

If  the  card  is  multiply  punched,  the  program  will  calculate  in  turn  the 
estimates  of  pollutant  concentrations  corresponding  to  the  specified  species. 

Card  3 will  provide  the  necessary  meteorological  data.  Of  the  entries 
made  on  this  card,  the  atmospheric  stability  index  (JSTAB)  should  be  deduced 
on  the  basis  of  wind  speed,  cloud  cover,  time  of  the  day  and  other  isolation 
parameters  according  to  the  scheme  outlined  by  Turner*-*.  Then  lid  height 
(NLID) , if  not  obtainable  by  processing  local  sounding  data,  is  estimated  on 
the  basis  of  JSTAB  and  the  monthly  average  afternoon  maximum  mixing  depths 
(meters  above  surface)  tabulated  by  Holzworth**.  The  twelve  (12)  average 
values  that  are  directly  coded  into  the  program  (FUNCTION  HMIX)  apply  to  the 
Washington,  D.C.,  region  and  are  ordered  by  month  (NMONTH)  as  follows:  400, 
570,  1000,  930,  1120,  1310,  1180,  990,  980,  570,  680,  and  480  (meters). 

Card  4 provides  an  accounting  of  the  receptors  (points  at  which  the 
model  must  estimate  the  pollutant  concentration),  the  aircraft  and  aircraft- 
engine  types,  and  the  airport-access  vehicle  roadways.  The  aircraft  area 
source  comprises  four  emission  sources:  (1)  ground  service  vehicles,  (2) 
auxiliary  poweruunits,  (3)  aircraft  taxi  within  the  area,  and  (4)  aircraft- 
engine  idle  within  the  area.  Therefore,  an  entry  should  exclude  the  aircraft 
area  source  (NACA"0)  whenever  the  effects  of  all  the  above  emission  sources 
are  suppressed.  The  non-aircraft  area  source  is  assigned  pollutant  emission 
levels  later  in  the  input  stream,  hut  mav  be  excluded  by  setting  NANA  to 
zero  (0). 


i 
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Card  5 provides  the  Indicators  for  user  selected  Input.  All  of  the 
punched  card  columns  require  that  all  corresponding  card  numbers,  punched 
with  all  associated  data,  appear  in  the  appropriate  order  of  the  input 
stream.  For  example,  columns  18  and  21  are  punched  in  order  to  indicate 
that  card  number  10  and  card  number  11  will  appear  in  the  input  stream 
(see  Table  3.1).  Both  of  theae  bears  entries  for  the  number  of  engines 
with  which  each  aircraft  type  is  equipped,  and  the  second  bears  entries 
that  Indicate  the  engine  types. 

Card  6 provides  the  location  of  specific  receptor  points. 

The  remainder  of  the  input  data  deck  is  described  in  a list  form 
in  order  to  identify  clearly  those  cards  that  are  required  for  every  appli- 
cation of  the  model  and  those  cards  that  are  used  to  over-ride  program 
default  values. 


Card 

7 thru 

Card 

11 

Card 

12 

Card 

13 

Card 

14 

Card 

15 

Ca  rd 

16 

thru 

Card 

18 

Card 

19 

Card 

20 

and 

Card 

21 

Card 

22 

Card 

23 

and 

24 

Card 

25 

thru 

i Card 

27 

Card 

28 

Card 

29 

Card 

30 

and 

Card 

31 

Card 

32 

Card 

33 

Ca  i u 

34 

Card 

35 

and 

Card 

36 

Card 

37 

Card 

38 

Refer  to  default  values  on  page  27. 

Refer  to  default  values  on  page  28. 

Refer  to  default  values  on  page  30. 

Airport  runway  coordinates  are  always  present. 

Refer  to  default  values  on  page  30. 

Refer  to  default  values  on  page  29. 

Refer  to  default  values  on  page  30. 

Airport  taxiway  coordinates  are  always  present. 

Refer  to  default  values  on  page  30. 

Runway  apron  coordinates  are  always  present. 

Refer  to  default  values  on  page  30. 

Is  required  for  access  vehicle  roadway  coordinates. 

Is  required  for  airport  terminal  area  coordinates. 
Refer  to  default  values  on  page  30. 

Refer  to  default  values  on  page  31. 

Refer  to  default  values  on  page  29. 

Refer  to  default  values  on  page  31. 

Refer  to  default  values  on  page  29. 

Is  required  for  airport  non-aircraft  area  coordinates. 
Is  required  for  airport  non-aircraft  area  pollutant 
emission  level. 


4 


Note  that,  for  the  hour  of  interest,  the  user  specifies  the  number 
of  arriving  aircraft  of  each  type  in  Card  7.  The  number  of  aircraft 
departing  during  that  same  hour  is  assumed  equal  to  the  number  of  arrivals. 
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fable  1.1.  Card  Input  Sequence 


Card  Type  Columns  Format 


Comment 


1 (l  card)*  1-80  (20A4) 

TITLE 


(1  card)* 

CO 

1-2 

(A4) 

THC 

5-7 

(A4) 

NOx 

9-11 

(A4) 

PART 

13-16 

(A4) 

(1  card)* 

HTAF.R0 

1-8 

F8.0 

WSP 

9-16 

F 8 . 0 

DIR 

17-24 

F8.0 

TEMP 

25-32 

F8.0 

3STAB 

33-40 

18 

SLID 

41-48 

18 

MONTH 

49-50 

12 

NR 

51-52 

12 

ll  card)* 

15 

NSR 

16-20 

T5 

NACT 

21-25 

15 

NENT 

26-30 

15 

NAVR 

31-35 

15 

NACA 

36-40 

15 

NANA 

41-45 

15 

(l  card)+ 

DFAL'LT  - 

1 9 

T 

2 

12 

3 

15 

4 

18 

5 

21 

6 

24 

7 

27 

8 

30 

9 

33  * 

10 

35-36 

11 

38-39 

12 

41-42 

13 

44-45 

14 

47-48 

15 

50-51 

16 

5 3-54 

Title  card  may  be  punched  or  left  blank. 


Four  fields  having  four  columns  each  are  provided 
so  that  one  or  more  pollutant  species  may  be 
chosen  for  the  model  rim. 


Height  (feet)  at  which  wind  speed  and  direction 
are  measured. 

Wind  speed  (knots). 

Wind  direction  (degrees). 

Ambient  temperature  (degrees  F). 

Atmospheric  stability  index. 

Lid  height  (tens  of  feet). 

Month  of  year  (required  whenever  NLID“0) . 

Hour  of  day  (required  whenever  NLID-0  and  JSTAB=4) . 

The  total  number  (£20)  of  receptor  locations  at 
which  the  model  must  estimate  pollutant  levels 
The  number  (£10)  of  different  aircraft  types 
(B727,  Ttf,  etc.)  using  airport  facilities. 

The  number  (<5)  of  different  aircraft  engine  types 
(JT8D,TPE331)  with  which  the  above  aircraft  are 
equipped. 

The  number  (<60)  of  motor  vehicle  roadways  leading 
into  the  airport. 

1 if  aircraft  area  source  is  included.  0 otherwise. 

1 if  non-aircraft  area  source  is  included.  0 otherwise. 

If  a blank  or  zero  appears  in  any  column  the 


default 

value 

is 

used 

• 

Card 

7 

data 

— 

user 

defined. 

Card 

8 

data 

— 

user 

defined. 

Card 

9 

data 

— 

user 

defined. 

Card 

10 

data 

— 

user 

defined. 

Card 

11 

data 

— 

user 

defined . 

Card 

12 

data 

— 

user 

defined. 

Card 

13 

data 

— 

user 

defined . 

Card 

15 

data 

— 

user 

defined. 

Card 

16 

data 

— 

user 

def ined . 

Card 

17 

data 

— 

user 

defined. 

Card 

18 

data 

— 

user 

defined . 

Card 

19 

data 

— 

user 

defined. 

Card 

22 

data 

— 

user 

defined. 

Card 

25 

data 

— 

user 

defined. 

Card 

26 

data 

— 

user 

defined. 

Card 

27 

data 

— 

user 

defined. 

If 


Table  3.1.  Card  Input  Sequence  (contd.) 


!ard  Type  Columns  Format 


Comment 


17 

56-57 

18 

59-60 

19 

62-63 

20 

65-66 

21 

68-69 

22 

71-72 

23 

74-75 

- 

- 

i (NSR  cards)* 

NRUSED 

1-2 

12 

XRECP 

3-10 

F8.0 

YRECP 

11-18 

F8.0 

ZRECP 

17-26 

F8.0 

7 (1  card) 

NAC(K) 

6-55 

(5X, 

, 1015) 

8 

(1  card) 

FLNDG(K) 

6-55 

(5X, 

10F5.0) 

9 

(1  card) 

FTKOF(K) 

6-55 

(5X, 

10F5.0) 

0 

(1  card) 

NGIN(K) 

6-55 

(5X.10I5) 

1 (1  card) 

INGN(K) 

6-55 

(5X.10I5) 

2 (NENT  x 4 cards) 

EMI  ( I , J ,K)  16-55 

(15X, 

5F10.0) 


Card  30  data  — user  defined. 
Card  31  data  — user  defined. 
Card  32  data  — user  defined. 
Card  33  data  — user  defined. 
Card  34  data  — user  defined. 
Card  35  data  — user  defined. 
Card  36  data  — user  defined. 


Any  reference  number. 

X coordinate  of  receptor  (mi 
Y coordinate  of  receptor  (mi 
Height  of  receptor  (ft 

NACT  entries  of  the  number  of  arriving 
aircraft  of  each  type  K.  (Number  of  de- 
partures assumed  to  be  the  same  as  the 
number  of  arrivals). 

NACT  entries  of  time  (hours)  in  landing 
mode  for  aircraft  type  K. 

NACT  entries  of  time  (hours)  in  take-off 
mode  for  aircraft  type  K. 

NACT  entries  for  the  number  of  engines 
with  which  aircraft  type  K is  equipped. 


NACT  entries  for  the  engine  type  with 
which  aircraft  type  K is  equipped. 


For  each  of  the  above  engine  types, 
enter  an  emission  rate  (lbs/hr)  for  each 
pollutant  K (specified  on  card  type  2). 
Cards  1 thru  4 for  each  engine  type 
throttle  setting  will  correspond  to 
taxi,  idle,  landing  and  take-off  modes, 
respectively. 


3 (1  card) 
DSRW 
HRW 


16-25  F10.0  Width  (mi)  of  initial  dispersion  on  runway. 

26-35  F10.0  Height  (mi)  of  initial  dispersion  on  runway. 
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Table  3.1.  Card  Input  Sequence  (contd.) 


Card  Type 

Columns 

Format 

Comment 

17 

56-57 

Card  30  data  — user  defined. 

18 

59-60 

Card  31  data  — user  defined. 

19 

62-63 

Card  32  data  — user  defined. 

20 

65-66 

Card  33  data  — user  defined. 

21 

68-69 

Card  34  data  — user  defined. 

22 

71-72 

Card  35  data  — user  defined. 

23 

74-75 

-1 

L 

Card  36  data  — user  defined. 

6 

(NSR  cards)* 
NRUSED 

1-2 

12 

Any  reference  number. 

XRECP 

3-10 

F8.0 

X coordinate  of  receptor  (mi 

YRECP 

11-18 

F8.0 

Y coordinate  of  receptor  (mi 

ZRECP 

17-26 

F8.0 

Height  of  receptor  (ft 

7 

(1  card) 
NAC(K) 

6-55 

CSX, 

, 1015) 

NACT  entries  of  the  number  of  arriving 
aircraft  of  each  type  K.  (Number  of  de- 
partures assumed  to  be  the  same  as  the 
number  of  arrivals) . 

8 

(1  card) 

FLNDG(K) 

6-55 

(5X, 

1 

NACT  entries  of  time  (hours)  in  landing 

10F5.0) 

mode  for  aircraft  type  K. 

9 

(1  card) 
FTKOF(K) 

6-55 

(5X , 

» 

NACT  entries  of  time  (hours)  in  take-off 

10F5.0) 

mode  for  aircraft  type  K. 

10 

(1  card) 

NGIN(K) 

6-55 

(5X.10I5) 

NACT  entries  for  the  number  of  engines 

with  which  aircraft  type  K is  equipped. 

11 

(1  card) 

INGN(K) 

6-55 

(5X.10I5) 

NACT  entries  for  the  engine  type  with 

which  aircraft  type  K is  equipped. 

12 

(NENT  x 4 cards) 

EMI  (I,J,K) 

16-55 

(15X, 

For  each  of  the  above  engine  types, 

5F10.0) 


enter  an  emission  rate  (lbs/hr)  for  each 
pollutant  K (specified  on  card  type  2). 
Cards  1 thru  4 for  each  engine  type 
throttle  setting  will  correspond  to 
taxi,  idle,  landing  and  take-off  modes, 
respectively. 


13  (1  card) 
DSRW 
HRW 


16-25 

26-35 


F10.0 

F10.0 


Width  (mi)  of  initial  dispersion  on  runway. 
Height  (mi)  of  initial  dispersion  on  runway. 


L 


I 

-'I 
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Table  3.1.  Card  Input  Sequence  (contd.) 


Card  Type 

Co lurans 

Format 

14  (1  card)* 

XI 

16-25 

F10.Q 

Y 1 

26-35 

F10.0 

Zl 

36-45 

F10.0 

X2 

*46-55 

F10.0 

Y2 

56-65 

F10. 0 

Z2 

66-75 

F10.0 

15  (l  card) 

VA1 

16-25 

F10.0 

VA2 

26-35 

F10 . 0 

VD1 

36-45 

F10.0 

VD2 

46-55 

F10.0 

TIME 

56-65 

F10.0 

TAIL 

66-75 

¥10.0 

lb  (1  card) 

FTAXI(K) 

6-55 

(5X, 

10F5.0) 

17 

( i card) 
FTXII(K) 

6-55 

(5X, 

10F5.0) 

18 

(1  card) 
FTXIO(K) 

6-55 

(5X, 

10F5.0) 

19 

(1  card) 
DSTW 

16-25 

F10.1 

HTW 

26-35 

F10.1 

20 

(1  card)* 

XI 

16-25 

F10.0 

Yl 

26-35 

Zl 

36-45 

X2 

46-55 

Y2 

56-o5 

Z2 

66-75 

Comment 


X coordinate  of  runway  (mi)  (Aircraft 
touch-down  end) 

Y coordinate  of  runway  (mi)  (Aircraft 
touch-down  end) 

Z height  (ft)  of  runway.  (Aircraft 
touch-down  end) 

X2  coordinate  of  opposite  end  point  (mi) 

Y2  coordinate  of  opposite  end  point  (mi) 

Z2  height  of  opposite  end  point  (ft) 

Initial  velocity  of  arriving  aircraft  (mi/hr) 
Final  velocity  of  arriving  aircraft  (mi/hr) 
Initial  velocity  of  departing  aircraft  (mi/hr) 
Final  velocity  of  departing  aircraft  (mi/hr) 
Aircraft  take-off  roll-time  (hrs) 

Exhaust  tail  length  (mi) 

NACT  entries  of  aircraft  taxi  speed  (mi/hr) 
while  in  the  gate  area  for  each  aircraft 
type  K. 

NACT  entries  of  aircraft  inbound  taxi 
speed  (mi/hr) 

NACT  entries  of  aircraft  outbound 
taxi  speed  (mi/hr) 

Width  (mi)  of  initial  dispersion  on  taxiway. 
Height  (mi)  of  initial  dispersion  on  taxiwav. 


X coordinate  of  inbound  taxiway  (mi) 

Y coordinate  of  inbound  taxiway  (mi) 

Z height  of  inbound  taxiway  (ft) 

X2  coordinate  of  inbound  taxiway  (mi) 
(Airport  terminal  end) 

Y2  coordinate  of  inbound  taxiway  (mi) 
(Airport  terminal  end) 

Z2  height  of  inbound  taxiway  (ft)  (Airport 
terminal  end) 
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Table  3. 

1 • Card 

Input  Sequence,  (contd.) 

Card  Type 

Columns 

Format 

Cotanent 

21 

(1  card)* 
XI 

16-25 

F10.0 

X coordinate  of  outbound  taxiway  (mi) 

Y 1 

26-35 

F10.0 

Y coordinate  of  outbound  taxiway  (mi) 

Z1 

36-45 

F10.0 

Z height  of  outbound  taxiway  (ft) 

X2 

46-55 

F10.0 

X2  coordinate  of  outbound  taxiway  (mi) 
(apron  end) 

Y2 

56-65 

F10.0 

Y2  coordinate  of  outbound  taxiway  (mi) 
(apron  end) 

12 

66-75 

F10 . 0 

Z2  height  of  outbound  taxiway  (ft) 

(apron  end) 

22 

(1  card) 
DSRA 

16-25 

F10.0 

Width  (mi)  of  initial  dispersion  on  apron. 

HRA 

26-35 

F10.0 

Height  (mi)  of  initial  dispersion  on  apron. 

23 

(1  card)* 
XI 

16-25 

F10.0 

X coordinate  of  inbound  apron  (mi) 

Y 1 

26-35 

F10.0 

Y coordinate  of  inbound  apron  (mi) 

Z1 

36-45 

F10.0 

Z height  of  inbound  apron  (ft) 

X2 

46-55 

F10.0 

X2  coordinate  of  inbound  apron  (mi) 

(from  here  we  start  the  inbound  taxiway) 

Y2 

Y2  coordinate  of  inbound  apron  (mi) 

(from  here  we  start  the  inbound  taxiway) 

Z2 

12  height  of  inbound  taxiway  (ft) 

24 

(1  card)* 
XI 

16-25 

F10.0 

X coordinate  of  outbound  apron  (mi) 

Y 1 

26-35 

F10.0 

Y coordinate  of  outbound  apron  (mi) 

Z1 

36-45 

F10.0 

Z height  of  outbound  apron  (ft) 

X2 

46-55 

F10.0 

X2  coordinate  of  outbound  apron  (mi) 

Y2 

56-65 

F10 .0 

Y2  coordinate  of  outbound  apron  (mi) 

12 

66-75 

F10.0 

Z2  height  of  outbound  apron  (ft) 

25 

(1  card) 

Omit  if  NAVR 

■ 0. 

EFUH(J) 

16-65 

( 15X, 

For  each  of  the  pollutants  specified  on 

5F10.0) 

card  type  2,  enter  the  urban  automobile 
emission  factor  (gm/km)  based  on  25  mi/hr 
traffic  speed. 

26 

(1  card) 

Omit  if  NAVR 

“0. 

EFUL(J) 

16-65 

(15X, 

For  each  of  the  pollutants  specified 

5F10.0) 

on  card  type  2,  enter  the  urban  automobile 
emission  factor  (gm/km)  based  on  10  mi/hr 

traffic  speed. 

27 

(1  card) 

Omit  if  NAVR 

- 0. 

DSAR 

16-25 

F10. 0 

Width  (mi)  of  initial  dispersion  on 
access  roadways. 

HAR 

26-35 

F10.0 

Height  (mi)  of  initial  dispersion  on 

L 


Card  Type  Columns  Format 


Comment 


28  (NAVR 

cards) * 

XI 

11-20 

F10.0 

Y1 

21-30 

F10.0 

21 

31-40 

F10.0 

X2 

41-50 

F10 .0 

Y2 

51-60 

F10.0 

Z2 

61-70 

F10.0 

VNOON 

71-78 

F8.0 

IFS 

79-80 

12 

29 

(1  card)* 

Omit  if  NACA 

= 0. 

XS 

16-25 

F10 . 0 

YS 

26-35 

F10.0 

STKH 

36-45 

F10.0 

WIT 

46-55 

F10.0 

30 

(1  card) 

Omit  if  NACA 

= 0. 

EFD(l) 

16-25 

F10 . 0 

EFD(2) 

26-35 

F10.0 

EFD(3) 

36-45 

F10.0 

EFD(4) 

46-55 

F10.0 

EFD(5) 

56-65 

F10.0 

31 

(1  card) 

Omit  if  NACA 

= 0. 

EFG(l) 

16-25 

F10.0 

EFG(2) 

26-35 

F10.0 

EFG(3) 

36-45 

F10.0 

EFG(4) 

46-55 

F10.0 

EFG(5) 

56-65 

F10.0 

32 

(14  cards) 

Omit  if  NACA 

= 0. 

SRVTIM(1,K) 

6-55 

33 

(1  card) 

Omit  if  NACA 

= 0. 

KAPU(K) 

6-55 

(5X, 

1015) 

34 

(1  card) 

Omit  if  NACA 

= 0. 

APU(J) 

16-65 

( 15X, 

5F10.0) 


Access  vehicle  roadways: 

X coordinate  of  first  end  point  (mi) 

Y coordinate  of  first  end  point  (mi) 

Z coordinate  of  first  end  point  (ft) 

X coordinate  of  second  end  point  (mi) 

Y coordinate  of  second  end  point  (mi) 
Z coordinate  of  second  end  point  (ft) 
Average  traffic  volume 

Roadway  classification:  l=congested, 
O=non-congested 


X coordinate  of  terminal  area  center  (mi) 
Y coordinate  of  terminal  area  center  (mi) 
Height  of  terminal  (ft) 

Side  length  of  terminal  (mi) 

Pollutant  emission  rate  of  diesel  engine 
powered  service  vehicles. 

CO  (gm/ gal) 

HC  (gm/gal) 

NOx  (gm/gal) 

PART  (gm/gal) 

SOx  (gm/gal) 

Pollutant  emission  rate  of  gasoline 
engine  powered  service  vehicles. 

CO  (gm/mi) 

HC  (gm/mi) 

NOx  (gm/mi) 

PART  (gm/mi) 

SOx  (gm/mi) 


NACT  entries  of  service  vehicle  operation 
time  (min)  during  the  aircraft  service 
operation  in  the  terminal  area.  See  Page  31 
for  a list  of  14  modeled  service  vehicles. 


NACT  entries  denoting  l=the  use  of  an 
auxiliary  power  unit:  0=no  auxiliary 
power  unit  used  for  aircraft  type  K. 


For  each  of  the  pollutants  specified  on 
card  type  2,  enter  the  emission  factors 
(lb/hr)  for  auxiliary  power  units. 


Tit  Mi’  3.1.  Cat'll  Input  Sequence  (contil.) 


Card  Type 

Column a 

Fo  rrna  t 

Comment 

»•>  ( l card) 

Omit  If  NACA 

- 0. 

F IDLE  00 

b-55 

(5X, 

NACT  entries  of  time  (hr)  for  aircraft 

10F5.0) 

K engine  idle. 

)«>  (1  card) 

Out  1 1 If  NACA 

- 0. 

TCND(K) 

b-55 

(5X, 

NACT  entries  of  time  (min)  for  aircraft 

10F5.0) 

K gate  occupancy. 

1/  (1  card)* 

Omit  If  NANA 

- 0. 

XS 

lb-25 

FI  0.0 

X coordinate  of  non-aircraft  area  source 
center  (ml) 

YS 

2b-)5 

F10.0 

Y coordinate  of  non-aircraft  area  source 
center  (ml) 

STKH 

tb-45 

FIO.  0 

Height  of  non-aircraft  area  source  (ft) 

WIT 

46-55 

FIO.O 

Side  tength  of  non-aircraft  area  source  (mi) 

IS  u card)* 

ihn t t If  NANA 

- 0. 

KMlT(.l) 

lb-b5 

(t«. 

For  each  of  the  pollutants  specified  on 

5F10.0) 

card  type  2,  enter  the  non-aircraft  area 
source  emission  rate  (lhs/hr) 

* 

Del  an  It  valu 

cs  for  these 

data  are 

not  available. 

See  Table  4. 

l.  Program  Default  Values  by  Card  Input  Sequence  Number . 

4 DESCRIPTION  OF  DEFAULT  DATA 


This  section  describes  for  the  user  what  is  required  as  input  data 
should  the  option  to  over-ride  program  constants  be  exericsed.  Also,  it 
lists  the  values  of  progmm  constants  (Table  4.1)  as  they  appear  in  the 
BLOCK  DATA  subprogram  of  the  computer  code.  All  or  some  of  the  values 
will  change  whenever  the  appropriate  card  is  inserted  into  the  input  stream. 

The  computer  program  is  designed  to  model  activity  of  ten  (10) 
different  aircraft  types.  These  are  listed  in  the  first  column  of  Table  4.1 
(page  27)  and  may  be  replaced  with  types  having  identical,  or  nearly  identical 
operational  parameters,  engine  configurations,  and  ground  service  requirements, 
ments. 

Operational  parameters  include  time  period  values  during  which  the 
aircraft-engine  throttle  setting  is  adjusted  to  one  of  the  following 
operational  modes: 

1 Landing  - Aircraft  touch-down  to  beginning  of  taxi  on  the 

inbound  apron. 

2 Take-off  - After  alignment  of  aircraft  with  runway  to  liftoff. 

3 Idle  - Arriving  aircraft  awaiting  gate  position. 

These  operating  times  are  shown  in  Table  4.1,  page  27,  columns  3 and  4, 
and  page  29,  column  6,  respectively.  Other  operational  parameters  presented 
on  page  29,  columns  2,  3 and  4 consider  the  speed  of  aircraft  while  they  are 
in  the  gate  area,  on  the  inbound  taxiway  and  on  the  outbound  taxiway, 
respectively.  Those  that  detail  the  average  performance  characteristics 
of  all  aircraft  types  during  flight  modes  of  operation  are  shown  on  page  30 
(Card  15) . 

Engine  configuration  and  their  emission  characteristics  are 

assigned  integer  values  to  denote  for  each  aircraft  type, and  each  engine 
type  the  number  of  engines  per  aircraft  and  their  pollutant  emission  rates 
during  4 modes  of  operation  (Taxi,  Idle,  Landing  and  Take-off).  These 
values  are  shown  in  Table  4.1,  page  27,  columns  5 and  b respectively,  with 
corresponding  pollutant  emission  rate  shown  on  page  28. 

The  ground  serivce  requirements  of  each  aircraft  type  are  fulfilled 
using  fourteen  (14)  different  pieces  of  motorized  equipment  (page  31,  Card  32), 
.ill  of  which  operate  within  the  gate  area  to  load  and  unload  cargo  and 
otherwise  prepare  the  airplane  for  its  next  departure.  It  can  be  seen 
l>y  the  service  times  entered  into  the  table  that  the  equipment  in  use  is 
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dependent  upon  the  type  of  aircraft  being  serviced.  The  pollutant  emission 
rates  (page  30,  Cards  30,  31)  are  presented  for  both  diesel  and  gasoline 
engine  powered  equipment.  Note  that  an  auxiliary  power  unit  will  provide 
electrical  power,  whenever  used  (Page  29,  Column  5),  for  the  entire  gate 
time  (Page  29,  Column  7)  of  the  aircraft.  Its  emissions  rates  are  shown 
on  Page  31,  Card  34. 


Table  4.1.  Program  Default  Values  by  Card  Input  Sequence  (Contd) 

Card  12  Pollutant  Emission  Rate  for  Each  Engine 
Type  During  4 Modes  of  Operation  (lbs/hr) 
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Card  13  Initial  dimensions  of  dispersion  on  runway. 

Width  » DSRW  - 0.030  miles 
Height  “ HRW  = 0.002  miles 


Card  15  Runway  parameters  for  arrival  and  departure  aircraft. 

VA1  * Runway  - arrival  initial  velocity  « 145  (mi/hr) 

VA2  » Runway  - arrival  final  velocity  » 25  (mi/hr) 

VD1  = Runway  - departure  initial  velocity  » 0.0  (mi/hr) 
VD2  » Runway  - departure  final  velocity  - 180  (mi/hr) 

TIME  - Take-off  roll  time  * .0111  (hrs) 

TAIL  ■ Exhaust  tail  length  = .8523  (mi) 


Card  19  Initial  dimensions  of  dispersion  on  taxiway. 

Width  = DSTW  - 0.030  (mi) 

Height  = HTW  = 0.002  (mi) 


Card  22  Initial  dimensions  of  dispersion  on  runway  apron. 

Width  = DSRA  = 0.095  (mi) 

Height  = HRA  = 0.002  (mi) 


Ca rd  25  Urban  automobile  pollutant  emission  factors  (gm/km)  based  on 
25  (mi/hr)  traffic  speed. 

CO  = 32.36,  HC  = 4.75,  NO  =3.46,  PART  - 0.19,  SO  =0.11 

’ x x 

Card  26  Urban  automobile  pollutant  emission  factors  (gm/km)  based  on 

10  (mi/hr)  traffic  speed. 

CO  = 70.18,  HC  =«  8.62,  NO  = 2.86,  PART  - 0.19,  SO  = 0.11 
’ x x 

Card  27  Initial  dimensions  of  dispersion  on  access  roadway. 

Width  = DSAR  =•  0.0095  (mi) 

Height  * HAR  - 0.001  (mi) 


Card  30  Pollutant  emission  factor  for  diesel  engine  powered  service 
vehicles  in  (gm/gal) . 

CO  = 126.6,  HC  = 21.9,  NO  - 185.82,  PART  = 5.9,  SO  - 0.0 

x x 


Card  31  Pollutant  emission  factor  for  gasoline  engine  powered  service 
vehicles  in  (gm/mi). 

CO  « 138.81,  HC  - 21.35,  NO  = 9.32,  PART  = 0.85,  SO  - 0.0 

X x 
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Table  4.1  (Cont'd) 


Card  16  Card  17  Card  18  Card  33  Card  35  Card  36 
Speed  Speed  Speed 

(ml/hr)  (ml/hr)  (ml/hr)  Time  (hr)  Time  (min) 


Aircraft  Type 

Gate 

area  taxi 

Inbound 

taxi 

Outbound 

taxi 

APU 

use  flags 

Engine 

idle 

Gate 

occupancy 

1 (Boeing  727) 

10 

15 

12 

1 

.033 

52. 

2 (Douglas  DC9) 

10 

15 

12 

1 

.033 

52. 

3 (Boeing  737) 

10 

15 

12 

1 

.033 

52. 

4 (Conviar  580) 

10 

15 

12 

1 

.033 

52. 

5 (BAC  111) 

10 

15 

12 

0 

.033 

52. 

6 (NAMCO  YS  11) 

10 

15 

12 

0 

.033 

52. 

7 (Beech  99) 

10 

15 

12 

0 

.033 

52. 

8 (Fairchild  FH  227) 

10 

15 

12 

0 

.033 

52. 

9 (Twin  Otter) 

10 

15 

12 

0 

.033 

52. 

10  (Piston  Engine) 

10 

15 

12 

0 

.033 

52. 
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Table  4.1  (Cont'd) 


Card  34  Pollutant  emission  factors  for  auxiliary  power  units  in  (lbs/hr). 
CO  - 2.82,  HC  - 0.11,  NO  - 1.24,  PART  - 0.0,  SO  - 0.0 

^ X 


Card  32 . Minutes  of  Service  Vehicle  Operation  While  Servicing 

Aircraft  Type  I 


Aircraft  Type 


Vehicle  Type 

1* 

727 

2* 
DC  9 

3* 

737 

4* 

C5 

5 

BAC 

6 

YS 

7 

B9 

8 

FH 

9 

TO 

10 

GA 

I Tractor 

66 

48 

85 

55 

50 

50 

0 

0 

0 

0 

2 Belt  Loader 

28 

15 

30 

0 

25 

25 

0 

0 

0 

0 

3 Container  Loader 

6 • 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 Cabin  Service 

12 

0 

15 

0 

0 

0 

0 

0 

0 

0 

5 Lavatory  Truck 

15 

15 

15 

10 

10 

10 

5 

5 

5 

0 

6 Water  Truck 

0 

10 

0 

10 

10 

10 

5 

5 

5 

0 

7 Food  Truck 

17 

17 

20 

10 

10 

10 

0 

0 

0 

0 

8 Fuel  Truck 

20 

15 

15 

10 

20 

20 

10 

10 

10 

0 

9 Tow  Tractor 

10 

5 

5 

5 

5 

5 

0 

0 

0 

0 

10  Conditioner 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

11  Airstart 

Transporting  and 

Diesel  Engines 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

12  Ground  Power 

Transporting  and 

Gasoline  Engines 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

13  Ground  Power  Unit 

Diesel  Engine 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

14  Transporter 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Also  serviced  by  an  Auxiliary  Power  Unit  (APU) 
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5 EXAMPLE  PROBLEM 


5.1  INTRODUCTION 

In  order  to  clarify  the  procedure  for  using  the  AVAP  Model  Abbreviated 
Version,  a "first-guess"  estimate  of  pollutant  concentration  is  presented 
for  the  example  airport-layout  shown  In  Figure  5.1  . it  is  instructive 
to  state  that  the  locations  of  line  type  sources  (runway,  taxiway,  apron, 
access  roadway)  are  specified  by  the  coordinates  at  their  edge  (from  end 
to  end),  but  the  locations  of  airport  terminal  and  non-aircraft  area 
sources  are  specified  by  the  coordinates  at  their  center.  Since  aircraft 
movement  on  the  runway  is  almost  always  into  the  wind  (given  as 
180  degrees)  the  pattern  of  runway-apron- taxiway  usage  emerges.  The 
pattern,  for  the  airport-layout  shown,  is  clockwise  starting  from  the 
northern  most  point  of  the  runway  (landing  aircraft  touch-down) . The 
input  coordinates  must  be  ordered  to  preserve  this  pattern.  Therefore, 
all  aircraft  prepared  for  departure  are  queued  before  take-off  on  Apron  2 
(outbound) . 

In  this  example,  of  the  ten  different  aircraft  types  for  which  the 
Abbreviated  Version  is  internally  coded,  seven  are  considered.  They 
are  listed  with  their  associated  engine  type  as  follows: 

(Boeing  727  & 737,  Douglas  DC9)  JT8D  turbofan  engine 

(Beech  99)  TPE  331  turboprop  engine 

(Convair  580,  NAMC  YS11)  All.  501-D13  turboprop  engine 

(Small  Piston  Engine  aircraft)  (>320  small  piston  engine 


5.2  STATEMENT  OF  THE  PROBLEM 

GIVEN:  Two  pollutant  species:  Carbon  Monoxide  and  Oxides  of  Nitrogen 


Height  of  aerovane 
Wind  speed 
Wind  direction 
Ambient  temperature 
Stability  index 
Lid  height 


10  ft 
13  knots 
180  degrees 
36  degrees  F 
4 

0 (Uses  coded  table  from  Holzworth) 


A 
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Receptor  locations 
Aircraft  types 

Engine  types 


Motor  vehicle  roadways 
Aircraft  area  source 
Non-aircraft  area  source 
Runway  end-point  coordinates  in  miles. 
Taxiway  end-point  coordinates  in  miles. 
Apron  end-point  coordinates  in  miles. 


3 (Coordinates  in  miles) 

10  (Use  defaults  provided  for 
these  types) 

4 (Use  defaults  provided  for 
these  types) 

5 (End-point  coordinates  in  miles) 
1 (Center  coordinates  in  miles) 

1 (Center  coordinates  in  miles) 


FIND:  The  "first-guess"  estimate  of  pollutant  concentration  at  receptors 

1,  2,  and  3 due  to  the  following  three  source  categories:  « 

1.  Aircraft  Area  Sources 

t 

Service  vehicles  used  in  the  terminal  area. 

Auxiliary  Power  Units  used  in  the  terminal  area. 

Aircraft  taxi  while  in  the  terminal  area. 

Aircraft-engine  idle  while  in  the  terminal  area. 

2.  Non-Aircraft  Area-Line  Sources 

Trial  area  source 
Access  vehicle  roadways 

3.  Aircraft  Line  Sources 

Runway  (Landing  & Take-off) 

Taxiways  (Inbound  & Outbound) 

Aprons  (Inbound,  Outbound  Queuing) 

Note  that  the  estimates  are  based  on  aircraft  arrival  counts 
entered  on  Card-type  7,  roadway  traffic  counts  entered  on  Card- 
type  28,  and  trial  area  source  emissions  entered  on  Card-type  38. 


5.3  SUMMARY  OF  THE  PROBLEM 

The  AVAP  Model  Abbreviated  Version  requires  at  least  twelve  input 
cards.  Depending  upon  the  number  of  receptor  points,  number  of  access 
vehicle  roadways  and  choice  of  default  values,  there  may  be  more.  The 
format  for  each  card  is  given  in  Table  3.1.  Table  5.1  lists  the  input 
for  the  sample  problem;  Table  5.2  lists  the  results.  Note  that  a descrip- 
tive list  of  program  values  starts  the  table  of  results.  Also  note 


if 


A 
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that  a tabulation  of  pollutant  emission  levels  for  each  area  and  line  type 
source  la  provided.  The  micrograms-per-cubic-meter  notation,  as  seen  in 
the  concentrations  tabulation,  is  used  for  Total  Suspendea  Particulates 
(FART)  if  it  is  one  of  the  pollutants  being  modeled. 


Table  5.1  Input  for  Example  Problem 
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Table  5.2  Results  of  Example  Problem 
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Table  5.2  Results  of  Example  Problem  (contd.) 
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Table  5.2  Results  of  Example  Problem  (contd.) 
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Table  5.2  Results  of  Example  Problem  (contd.) 
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APPENDIX  A 

CONCENTRATION  CONVERSIONS 
■j  4 p 

Concentration  (ug/m  ) = 1.22  x 10  ^ • MW  • concentration  (ppm)  where 

P = pressure  in  atmospheres, 

T = absolute  temperature  (°K) 

MW  = molecular  weight 

Example:  N02 

P = 1 atm 

T = 77°F  = 25°C  = 273  + 25  = 298°K 
MW  = 14  + 32  = 46 

Concentration  (yg/m3)  = 1.22  x 104  x x 46  x concentration  (ppm) 

= 1883  x concentration  (ppm) 


for  T = 32°F, 

7 

Concentration  (ug/m  ) = 2056  x concentration  (ppm) 


Note  that  for  a pollutant  class  like  THC  or  NO  , it  is  customary 

7 * 

to  convert  from  ppm  to  yg/m  by  representing  the  class  in  terns  of  a single 
pollutant  such  as  CH^  or  N02>  Consequently,  under  the  same  conditions  as 
given  above  (T  = 32°F) , 

7 

THC  Concentration  (yg/m  ) (treated  as  CH^)  = 


1.22  x 10^  x 1 x 16 

m 


715  x TOC  cone,  (ppm) 


and 


N0x  Concentration  (yg/m3)  (treated  as  N02)  = 2056  x NO^  cone,  (ppm) 


START 


APPENDIX  B 


FLOW  DIAGRAMS 


READ  TITLE 

READ  POLLUTANT  CHOICE 


SCAN  CARD  TYPE  2 FOR  PROPER 
SPECIFICATION  OF  POLLUTANT (S) 


NPT  = Number  of  pollutants 
chosen. 


READ  CARD  TYPE  3 


HTAERO,  WSP,  DIR,  TEMP 
JSTAB.  NLtD 


READ  CARD  TYPE  4 

HSR,  NACT,  NENT,  NAVR, 
NACA,  NANA 


READ  CARD  TYPE  5 
DFAULT  (N),  N=1 , 23 


READ  CARD  TYPE  6 
R,  XRECPr,  YRECPr,  ZRECPr 


O 


B-fc 


B-7 


READ  CARD  TYPE  22 
DSRA,  HRA 


DFAULT (13) :0 


DSRA  = 

.095 

HTW  = 

002 

r 


COMPUTE  POLLUTANT  J EMISSION  DUE  TO  AIRCRAFT 
IDLE 

EMIDLj  = lNACK  x NGINk  x EMI  j lj  x FIDLEk 

K 

where  L = 2 


COMPUTE  POLLUTANT  J EMISSION  DUE  TO  AIRCRAFT 
GATE-AREA  TAXI. 

TERHJ  * \ » nack  x ngink  X emi,  LJ 


COMPUTE  POLLUTANT  J EMISSION  DUE  TO  AUXILIARY 
POWER  UNIT  USE 

EAPU , = APU,  I NAC,  x TGNDK  x KAPU„ 

J J K K “607"  K 


TOTAL  EMISSION  OF  POLLUTANT  J IN  THE 
TERMINAL  AREA 

EMITj  = PLLTNTj  + EMIDLj  + TERMj  + EAPUj 


NANA:  0 


^ READ  CARD  TYPE  37 
X,  Y,  STKH,  WIT 


READ  CARD  TYPE  38 


KSTAB 

= 0 

HLID^>i- 

[< 

KSTAB  = 1 


DLZ  = 7 

t 

DSS  = ^ 
KL  = 0 


2.  x STKH 
2.4  x 5280. 


SINWD  = SIN  (WDRAD) 
COSWD  = COS  (WDRAD) 

CWS  = X' COSWD  - Y- SINWD 

DWS  = X-SINWD  - Y- COSWD 

CWR  = XRECP- COSWD 
-YRECP-SINWD 

DWR  = -YRECP-SINWD 
-YRECP -COSWD 


COMPUTE  CROSS  WIND  AND  BOWNWIND  DISTANCES  FOR 
SOURCE-RECEPTOR  PAIR 

CW  = CWR  - CWS 

OW  = DWR  - DWS 


^/'THE  RECE^\ 
TOR  UPWIND  OF  THE 
sSSOllRCE . DW 
\WIT.>^ 


PFA  = 1. 

I BACK  = 0 
NRFLAG  - 0 


K 14 


/IS  TMCS. 
/RECEPTOR  \ 
/UPWIND  OF  THE  \ 

Downwind  edge  of  the 

.SOURCE? 


/ns  THE  \ > 

/CRITICAL  LENGTH\. 
^UPWIND  FROM  THE  RE.CEP-^^ 
TOR  DOWNWIND  OF  THE  DOWNWIND 
EDGE  OF  THE  SOURCE? 

v 'WR-XL)  (DWS+WIT)  / 


/niOES  THE  \ 
/CRITICAL  LENGTH 
^UPWIND  FROM  THE  RECEP-  ^ 
TOR  EXTEND  BEYOND  THE  UPWIND 
EDGE  OF  THE  SOURCE? 

'V>SE[)WR-XL)  < (DWS-W1T ) Z' 


l BACK  - 1 


m 
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FUNCTION  CAVL 

Purpose: 

To  compute  the  coupling  coefficient  due  to  a line  source  of  arbitrary 
spatial  orientation,  at  the  receptor. 


-End-point  coordinates  of  line 
source  (mi)  (X-axis  chosen  to 
be  along  wind  vector) . 

-Receptor  coordinates  (mi). 

-Wind  speed  (mi/hr);  Stability; 

Mixing  height  (ft);  Critical 
distance  for  vertical  mixing 

(mi) ; Pseudo  downwind  distance  \f 

for  horizontal  spread  of  line 
source  (mi) ; Pseudo  downwind 
distance  for  vertical  spread 
of  line  source  (mi). 

-Constant  coefficients  used  in 
the  line  source  dispersion 
equations. 

Output: 

CAM  (XR,  YR,  2R) Coupling  coefficient  (mi  ^)  at 

the  receptor  point  (XR,  YR,  2R) . 

Procedure: 

1.  Test  whether  the  receptor  is  located  with  respect  to  the  line  source 
such  that  the  concentration  is  completely  negligible. 

2.  If  angle  between  wind  vector  and  line  is  sufficiently  small,  and 
line  source  is  sufficiently  long,  a flag  is  set  for  the  line  to  be 
segmented.  Each  segment  is  then  treated  as  an  individual  line. 

3.  Compute  effective  downwind  distance. 

4.  Compute  horizontal  and  vertical  dispersion  coefficients,  using 
external  function  routines  SIGY  and  SIGZ. 

5.  Compute  the  Z-component  of  the  dispersion  expression. 


Input: 

XW1,  YW1 , ZW1;  XW2,  ZW2 


XR,  YR,  ZR 

WS;  JSTAB;  ULID;  XZ;  SUDOY;  SUDOZ- 


C0EF1 ; C0EF2- 


1 
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b.  Test  whether  the  line  source  has  a uniform  density.  If  it  is  a 
runway  used  for  aircraft  arrival  or  departure  (nonuniform  line 
density),  subroutine  QM1D  is  called. 

7.  Compute  and  output  the  concentration  for  the  given  receptor. 

functions  Called: 

SI GY,  SIGZ 

Subroutine  Called: 

QMOD 


Special  Function  Called: 

CRT  (X) 

$ 


COMPUTE  DOWNWIND  DISTANCE: 

DWI>XR-X 
DWDY =DWD+SUDOY 
DWDZ=DWEH-SUDOZ 


CALL  SIGY,  SIGZ  TO  COMPUTE 
DISPERSION  COEFF.  SIGH,  SIGV 


COMPUTE  EXPRESSIONS  TO 
BE  USED  IN  THE  ARGUMENTS 
OF  ERROR  FUNCTION  FOR 
HORIZONTAL  DISPERSION: 
AL;  BA1;  BA2 


COMPUTE  ARGUMENTS  OF 
EXPONENTIAL  FUNCTIONS 
FOR  VERTICAL  DISPERSION: 
ARGZ1;  ARGZ2  


SIN<|»|  =0.1 
SIN6I :0.1 


OOMPUTE  ARGUMENTS  OF 
EXPONENTIAL  RINCTIONS  FOR 
VERTICAL  DISPERSION  USING 
SMALL  ANGLE  APPROXIMATION 


a 


B-22 


B-27 


FAC3=0. 5*SIGV/ 
(C0EF1*HLIEM) 


DWD: 2*XZ 


1 

< 

DENZL=0.47*HLIDM 

DENZL2=DENZL**2 

FAC  1= FAC  3 
FAC2=FAC3 


JSIN*#. I :0.1 
ISINOI : 0 . 1 


COMPUTE  Z -ARGUMENTS 
ARGZ1 , ARGZ2  USING 
34ALL-ANGLE 
APPROXIMATION 


CGMPUTE  Z-ARGUMENTS 
ARGZ1 , ARGZ2  USING 
GENERAL  FORMULA 


FAC1=EXP(ARGZ1) 

FAC2=EXP(ARGZ2) 

FAC1-FAC1+ (DWD- XZ) * (FAC.V  FAC1 ) /XZ 
FAC2bFAC2+ (EWD-XZ) * (FAC3 -FAC2) /XZ 


E 


4 
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SUBROUTINE  QMOD 


Purpose: 

To  compute  the  linear  distribution  (in  (Length)  *)  of  pollutant 
along  runway  due  to  aircraft  emission  during  landing  or  takeoff. 


Input : 

YS1 Distance  along  runway  measured  from  the  tip  of 

exhaust  plume  near  the  starting  end  of  runway. 

TAIL  Length  or  penetration  of  the  exhaust  plume  of 

aircraft  at  rest. 

DL Total  length  of  the  smoke  slug  on  the  runway. 

A Acceleration  (or  deceleration)  of  the  aircraft. 

V12  Initial  Velocity,  VI,  squared. 

VS  Average  velocity  of  the  exhaust  particles  relative 

to  the  air  mass  in  the  tail  or  exhaust  plume. 

WS2 Wind  speed  squared. 

WSC  2* (Wind  Speed)* (-cosine  of  angle  between  runway  and 

wind  vector) . 

RR  A/G,  where  A is  the  acceleration  of  the  aircraft  and 

G is  the  normalization  constant  for  line  density. 


Output : 

QL 


V 


Procedure : 

1.  Convert  the  quantity  YS1  to  the  distance  measured  from  the 
physical  end  of  runway,  i.e. , XS1. 

2.  Use  the  line  density  formula  to  compute  QL(1/Length) . 


i 


-30 
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SUBROUTINE  QM3D  (Continued) 


> 


B-33 


r 


FUNCTION  TRAN 


Purpose: 


To  compute  the  coupling  coefficient  at  the  receptor  point  due  to 
a point  or  area  source. 

Input: 

1.  Meteorological  parameters:  wind  speed;  stability;  mixing  height; 
critical  distance  for  mixing. 

2.  Source  parameters:  horizontal  and  vertical  spreads;  pseudo  trans- 
port times  corresponding  to  the  horizontal  and  vertical  spreads; 
area  source  flags:  KSTAB,  NRFLAG,  and  IBACK. 

3.  Receptor  parameters:  downwind  and  crosswind  distances;  receptor 
he ight . 

Output : 


Point  or  area  source  coupling  coefficient:  TRAN. 

Procedure: 

1.  If  the  effective  stack  height  exceeds  the  height  of  the  sky  lid, 
then  stability  index  is  reassigned  according  to  the  current 
hourly  atmospheric  stability  and  the  source  flag  KSTAB  conputed 
in  PLUME. 

2.  For  point  source  and  area  source  with  source  flag  NRFLAG  = 0, 
compute  the  travel  time  for  z dispersion  from  the  center  and 
that  for  y dispersion  from  the  downwind  edge  of  the  source. 

3.  For  sources  with  NRFLAG  = 0 the  effects  of  ground  and  sky  lid 
are  treated  by  the  image  method.  Up  to  6 terms  are  included  in 
the  coupling  coefficient. 

4.  For  area  source  with  NRFLAG  = 1,  the  travel  times  from  the  down- 
wind and  upwind  edges  of  the  area  sources  are  determined  on  the 
basis  of  receptor  location  relative  to  the  area  source.  These 
plus  the  pseudo  travel  time  1 due  to  the  Z- spread  are  used  to 
compute  the  Z-dispersion  coefficients  CT^ (T^ ) andC~z(T?). 

5.  For  area  source  with  NRFLAG  = 1,  the  y-dispersion  coefficient 
(Ty(TT)  is  determined  on  the  basis  of  the  pseudo  travel  time  Ty 
due  to  the  y- spread  and  the  travel  time  from  the  downwind  edge 
to  the  receptor. 

6.  The  coupling  coefficient  for  area  source  with  NRFLAG  = 1 is  then 
computed  using  the  integrated  expression  for  "near"  source. 


if 


I 
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7.  If  the  source  flag  IBACK  is  1,  part  of  the  area  source  is  to  be  treated 
as  "near"  and  part  as  "far"  area  sources.  When  both  contributions  to 
the  coupling  coefficient  are  computed  and  summed,  IBACK  is  then  set  to 
0. 

Functions  Called: 

SIGY,  SIG2 


[ 
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E3=EXP(EXPZ3) 
ADD3=ADD1*E3/E1 
Z4=Z3+2.*ZKEEP 
EXPZ4=- . 5* (Z4/SZ) 2 


E4=EXP(EXPZ4) 

ADD4=ADD1*E4/E1 


1 


ADDS=ADD3 

ADD6=ADD4 


E5  = EXP (EXPZ5) 

ADD5  = ADD1*E5/E1 
Z6  = Z5  + 2.*ZKEEP  ? 
EXPZ6  = - . 5*(Z6/SZj 


I BACK: 1 


CALCULATE  T2  AS  TRANSIT 

tine  from  upwind  edge 

OF  SOURCE: 

T2  * TZ  + (DW+WITVWSAVE 


IS  THE  \ 
RECEPTOR  DOIVN- 
WIND  OF  THE  DOWNWIND 
S0)GE  OF  Tin:  SOURCE?^ 
\fDW-WIT)>0? 


CALCULATE  T2  AS  TRANSIT 
TIME  OVER  THE  CRITICAL 
LENCTI 1,  XL: 

T2  = TZ  + XL/WSAVE 


/ IS  THE  \ 
/ receftor  down-^' 

WIND  OF  THI-:  DOWNWIND 
sJBDGE  OF  TOE  SOURCE.? 
XlDW+WIT-DELY)  ' 
\ >0? 


CALCULATE  T1  AS  TRANSIT 
TIME  FROM  DOWNWIND  EDGE 
OF  SOURCE  AND  SIMILARLY, 
FOR  TT  AS  TRANSIT  TIM: 
FROM  DOWNWIND  EDGE  FOR 
Y DISPERSION: 

T1  = TZ  + (DW-WIT)/WSAVE 
IT  * TY  + (DW-W7~)/WSAVE 


SHORTEST  TRANSIT 
TIME  FROM  ANY 
PORTION  OF  SOURCE: 
T1  = TZ 
^ _ ty 


DELZ/2-4 


SIG2  (JSTABT.Tl ) 
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APPENDIX  C 

FORTRAN  LISTING  OF  THE  AVAP  MODEL  ABBREVIATED  VERSION 

C...  AVAP  MAIM  FGM  OOOOOOOl 

R*A(*A  TITLE  1(20), UTLF(IO)  00000002 

H E AL  * A NAME ( 10)  00000003 

RFAL*b  SRNAMF (1A) ,NAUT(5 I OOOOOOOA 

INTEGER  OUTPUT, 0FAULT(23>  00000005 

DIMENSION  F M I T ( 5,2)  00000006 

D IMF  NS  I UN  I K A( 101 , I KB ( 10)  00000007 

DIMENSION  SD Y ( 2 I 00000006 

( c M ( T CONTAINS  FIVE  POLLUTANT  FMISSIEN  PATES  f-  OR  TWO  APM  SOUPCFS  00000009 

DIMFNSION  XRECP(20),YRECP(20),2RFCP(20),NPUSFC'I20)  OOOOOOIO 

DIMENSION  LCMI10)  00000011 

DIMENSION  X A( 7) ,YA( 7) ,ZA( 7) , XR  ( 71  ,YB( 7)  ,Z8(  7)  00000012 

DIMENSION  FTXY(2I ,EMLN(5) ,FMTF(5)  00000013 

DIMENSION  CCUPI  20)  ,TERM(5)  .EMTXK  5 ),EMTXfl<5  ),  0000001A 

- RwICLI 5,2) ,FMIO( 5) ,EAPU( 5)  00000015 

0 iMfcNSI ON  A L SC ( 60 , 6 ) OOOOOOlfc 

DIMENSION  A L SF ( 60 , 6 ) 00)00017 

DIMENSION  PI LTNT( 5) , SRVHF ( 1A ) 00000016 

DIMENSION  CWS(2 ) ,0WS(2) , C WR ( 20  I , DWP ( 20 1 , P FC DA T ( A, 5, 20)  00000015 

DIMENSION  E'EFFI  5)  .KSTABI5)  00000020 

u mmon/oeflt/nac( ioi .flndgi ioi , ftkdfi ioi.ngini ioi , ingni ioi , 00000021 

lPMI(5,A,51,rSRw,HRW,  FT  AX  I ( 10)  , 00000022 

2ETXII (10) ,FTXIL)( 101 ,OSTW,HTW,DSPA ,HRA,FF()H( 5) *EFUl (51 ,D SAP, HAP,  00000023 

3SPVTIM( 10,1A),KAPU( 10I,APU(5),FIDLE ( 10) , TGND( 10 ) , DF AULT  0000002  A 

COMMC N/OELTA/DFLY.DELZ  00000025 

CCMMDN/FUEI /FFD(5) , EFG( 5 ) 00000026 

CCMMDN/LN/Dl , XW1 , YW1 ,ZW1 .XW2.YW2, ZW2.CDFF 1 ,C0FF2, VA1 , VA2,VT l,V02,  00000027 
-C  ,T 1 ME  , V A 1 2 ,VA22, VD12,VD22,WS2,WSC, JAC , SNAN , C S AN, VI, V2, VI 2, V22,  COO 00028 

-T AIL , VS.RR , SP  00000029 

CCMMON/LM/  I ACC  00000030 

rf MMON/LGC/DW.CW  00000031 

f CMMt  IN/MET  / W S , WO , JST  AH , ML  l[),TFMP,Xl  , SUDC1Y , SUCPZ  00)0003  2 

C C 3MUN/PI 3/ X S ( 5 I ,YS(5» , ST  KH ( 5 ) , W I T ( 5 ) C000003  3 

C CMMON/PGL/NPT, IJ (5 ) 000)00?A 

f CMMCN/RECPT/HRECPT ,HTAFRn,ZRFCPG  G0000035 

f CMMDN/PI SE/ZSS  00000036 

C CMMGN/SE UO C/TXT, TYT.TZT  00000037 

TMMON/Sf RVHl /SONAME  D0000038 

( CMMON/XTPAN/WSMO,NCALM,SUTOPI  00)00039 

DATA  NAME /• CD  ' , ' T HC  *,'NOX  • , 'PART' ,• St  X ',5*'  '/  OOOOOOAO 

:ATA  NAUT/'Cf  (PPM)  • , *THC ( PPM  1 • , *N  1X(PPM 1 • , •PTIM/f M) • , ' SD2( PPM) •/  OOOOOOR 1 
DATA  BLANK/'  '/  0000004? 

^ A T A S02 P I /2. 5066283  / 0009004-* 

C CCNVF  CONVERTS  LHS/MI**3  INTO  MICRCJ  GM/M**3  OOOOOOAA 

C A) A CDNVF/ . 10882 139/  000000A5 

DATA  ClN/1250.,666. ,205A. ,1 • ,285A.,5*l./  OOOOOOA6 

INPUT  = 5 000  0 )0A  7 

DUTPUT  = 6 000000A8 

RFAU(  INPUT,  1 71  TITLE1  OOOOOOA9 

PPADI INPUT , 19)  (TITLFII ) . T-1,10)  C0000050 

WP1  IF  (OUTPUT,  17)  TITLE-1  00000051 

WPI TE (CUT PUT, 31  I ( T I T L E ( I ) , 1 = 1,  101  00000052 

17  F DR  M A T ( 20A A ) 00000053 
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19  ‘ f PM7..T  l 10A5  ) 

S'  AN  TK  POLLUTANT  SPECIFICATION  LARI'. 

NP1  = 0 
I = 1 
l = l 
If l AG  = l 

5 )0  IF  (NAPt  < l I.FQ.T!  HEI  L»)  CCj  TO  50 2 
SOS  l =1*1 

if  1 1 .u.  io)  r,c  to  500 
I * l+l 

IHI.oT.10)  GO  TO  SOI 
l = 1 

c.r  io  s oo 

5c  l in  If  lAG.FO.n  GO  TO  503 
c.r  T(j  50s 

S I F|  11  Tlf  (L  ) .rC.'iLANK)  GO  TO  505 
\FT  = MT»i 

l HAG  = 0 
505  IJ(NPl)  = I 

IF INPT.rw.l  I GO  T 0 510 
I n I Jlffl  | .!  C.I  J(NPT-l)  ) GO  T ('  507 
510  ( TNI  I Mi1- 

I FI  1 .1  T.  10)  r.n  T 3 505 
OC  TL  505 

C0J  v»F  I If  (OUTPUT  * 50a  1 ( TITIF1  I ) • I = 1 * 1 0 ) 

SO”  c-Ch,'«AT  ( 1 HI , I FRPPO POLLUTANT  NAPF  NOT  FOUND  * , 10A5  ) 

G(  TC  505 

50  7 »al TF(COl PUT ,509)  TITLIIL) 

505  frFVMIlHl,'  DU°l  ICATF  POLLUTANT  PFGUFST  ',A51 

\pi  — r.pi-i 
GC  T I 505 
‘15  CCMIM.'f 
f 1 a I = 7 

' 1 5 V f - 15 

r.c„Y  = i 
\ A Pfc  = 7 

(..  JfJi ..  F-FIl.IT  OF  AFFOVANE  FOR  WINO  MFASUPCMFNT 

P f Ai,(  INPUT ,5713 ) HT AFRO ,KSP, 01 P ,TEPP, JSTAH.NL IP.NMONTH, NP 
5713  f CPMA  T io  a.  0,21  8,21  71 

p F A l ( INPUT,  12  ) NSF,NACT,FFNT,NAVP,NACA,NANA 

..  p-IFCT  offault  options. 

• f Al  ( INPUT,  1 ) (PFAULTU  >,  1 = 1,231 

C 

r . . . EAT  F F f F PT  r P CIOPDINATFS  X,Y,  IN  MILES,  l IN  FEET. 

-‘•AU  ( If. PUT,  5 711  I ( NPUSFC  (Nl  .XPfCPIN)  , YP  ECP  I N ) , l RF  t P ( N ) ,N=  1 , NSP  ) 
Or  77C l a 1*1, 5 
•'C  22  63  II  = 1 , N°  T 
*7  = IJILL I 
C 22F  A a 3=1,  NSP 
225’  " CU C A II  a.  i ,k2, NT  1-0.0 
T r y - T | M O 

«r  - i ip 
•'*  - aSp»1.15 
J S I 5 H F - JS  TAh 


00000055 
00000055 
00000056 
00000057 
00000058 
00000059 
00000060 
00000061 
0000006 2 
00000063 
00000065 
00000065 
00000066 
00000067 
00000068 
00000069 
00000070 
00000071 
00000072 
00000073 
00000075 
00000075 
00000076 
00000077 
00000078 
00000079 
00000080 
ooooooei 
00000062 
00000083 
00000085 
00000085 
00300086 
00000087 
00000088 
00000089 
00000090 
00000091 
00000092 
00000093 
00000095 
00000095 
00000096 
00000097 
00000098 
00000099 
00000100 
00000101 
00000102 
COO 001 0 3 
00000105 
00000105 
00000106 
00000107 
00000108 


JHIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
IJKNiOUrY  EUaJUSHU)  XO  DDC 
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FUNCTION  HM|X  WILL  RETURN  a VALUE  FOR  MIXING  DF  P TH  00100109 

HI IL  = M 10*10  00)00110 

fflHL ir.FO.O.OI  HI !P=HMJX ( JSTABB,NMCNTH,NR,WS)  00000111 

WRI TF (OUTPUT , 34)  NSP * NR W Y , NT  X 1 ,NAPF  ,NAC A , NACT , NENT , NS VP  ,NAVR  .NANA  00000112 
WP1  TF  (OUTPUT  ,r>71 1 ) N$P  , ( XPEC P t N ) , YR FCP  ( N)  , ZRF  CP  ( N ) , N*  1 , NSR  I 00000113 

*P I T E (OUTPUT ,97 44  I (NPUSFP(N1,N*1,NSP>  00000114 

WPI Tf  (OUTPUT, 9734  1 HT AF PC, W S ,0  I R , T F MP , JST AB , HL 1 D 00000115 

COMPUTE  SINE  A NO  COSINE  OF  THE  Wl  NC-  ANGLE  00000116 

WCPAt =W0*. 01745  00000117 

S 1 Nwf  = S I MWPPAO  I 00000118 

LFSWC*f OS(WOPAO)  00000115 

A =wS*  S l NWO  CO  TOO  1 20 

B = W S*C ( SWD  00000121 

HL  = Hl  10*0. 4 7/5280.  00000122 

PS=0.  00)00123 

KL=2  00000124 

T X-0.  00000125 

TY  = 0.  00000126 

T Z = 0 . 00000127 

FALL  PSFUPniDS, A.H.TX.TV ,TZ,l)X,OY, JSTAb.HL , KL ) 00000128 

XL=WS*TZ  00000129 

W S2  = WS*WS  00000130 

« an  = 270. -WP  000  )0131 

W AN  = WAN*3 . 141 5927/ 180.  00000132 

LEAN  = CGSIWANI  00000133 

SNAM*  SIM  WAN)  00)001^4 

INPUT  AIRCRAFT  AWRIVAI  ANT  L)E  P APTURF  ACTIVITY  BY  TYPE  . 00000135 


IFJCFAUl T ( 1 ).NE.O) 

1RFADI INPUT, 5) (NAC( I ) , 1 = 1, NAC  T ) 

WHITE (OUTPUT, 5o0)  NACT, If  AC  ( K ) , K=  l , N ACT  I 
hi)  (CPMAT  ( IX,* ARRIVAL  ACTIVITY  FOP  AIPCRAFT  TYPES  1 THROUGH*, 13 
-.3X.10I3) 

INPUT  TlMFlUfUE)  oPCNT  If  LANDING  MOFE. 

fF(OFAULT (2  I .Nt.O) 

1 E F AT ( INPUT, 25)(FLNDG(  I), 1 = 1, NACT) 

„■  c 1 T F (OUTPUT,  59  7)  (FLfCG(I),I=l,NACT) 

INPUT  TIME! HOUR ) SPFNT  IN  TAKE-CCF  MCOE. 

I F ( l FAUITJ3I.NE.0) 

1 R F A I ( INPUT , 2 5) ( c TKOF ( I ) , 1*1, NACT) 

WRITE (OUTPUT,  59 81  ( F TKOF ( i I, 1 = 1, NACT  I 


00000136 
00000137 
OOOOOD8 
00000129 
00000140 
00000141 
C0000142 
00000143 
00000144 
0000014* 
00000146 
00)00147 
000001 4 E 
000)0145 
00  )0014  " 
000  101*1 
000)01 * 
OOCOOl*-4 
000)01 64 


! NF  o T TfF  NUMBER  IF  ENGINES  USED  WITH  FACH  AIRCRAFT  TYPE.  00)00165 


00000156 

!F(LFAULT(4).NF.O)  00000157 

1C FAI  I INPUT, 5) (NO  I N(  I ) , I =1 , NACT ) 000001 5F 

000)01*° 

INPUT  IK  FNG1NF  TYPE  I NT  F X OF  THE  ABCVF  ENGINFS.  00000160 

0000016 1 

I F (l)FAUL  T (S  » .Nfc  .0  ) 0000016  2 

IF  F AT  C IN'>UT,  5 M INGN(  I)  , I * 1 , NACT  I 00)00163 
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v • • 

input  tip  tut  sfecmet  pihiutanttsi,  the  emission  ratmir/hpi 

r..  UF  EACH  AIRCTAfT  -ENGINE  TYPt  FCIK  THE  FOLLOWING  FOUR  (4)  KfOFS 
IP  CPfPAT  ll'N : 1.  TAXI.  2 • IEIE.  3.  LANDING,  A.  TAKE-OPF. 

IF(TF  AUl  HtI.F3.OI  GC  TC)  556 

C .. 

DO  55b  I - 1 » Nc  NT 
PC  555  J--1.A 

55->  P F Ai  ( INPUT  ,3  ) <1-  MI  l I , J,  I J(  K)  ) ,K*  L,  NPT) 

*»56  CTNTIKIF 
IF  = 5 

C 

I FI  TP  Alii  T ( ) .FO.OI  GO  TO  55  7 

1 5 = NF  NT 
C 

,'W  'C  ‘.Ad  I = 1,1b 
PC  ‘♦At  J = 1 , A 

*P  I Tu  (C  UT  PUT  , lb  3)  ( J,  (MUI  I,  J,  I J(M»  ,K=1  ,NPT  I) 

A Ad  U C N T 1 NUF 

r . . 

i ..  INPUT  INITIAL  "TENSIONS  OF  RIJNWAY  WIDTH  AND  HEIGHT. 

* • • 

IF  (,)P  AUl  T 171  ,NE  .01 
I r AO  I !N°UT,  3 ? 3 1 OSRW.HFW 
- . St-  w - DFRW/^.A 
MR  „ = HP  w / 1 • ? 

RL  - 0 

TAU  F Sr  UP  C IPSFw.A.d.TXl.TYl.TZl.rX.DY.JSTAb.HPW.KL) 

Si UOY  = WS  * T X 1 
SUDl/  = *S  * TZ1 
i . . INPUT  P UP  « AY  CCDR31NATFS. 

r # 

° c 4U  I INPUT ,1255)  XI,Y1,Z1,X?«Y2,22 
wL  I IE ICUTPUT  ,971b)  X1,Y1,Z1,X2,Y2,Z2 
Z 1 = 21/5280. 

12  = id 526C. 

JINX  = X? 

-LNY  - Y ? 


I 3UT  EUNWAy  PAPANFTFRS:  I A I T I A t AND  FINAL  VF10CITIFS  DUPING 
,-FIVAl  AM  TEPAUTURF.,  T A K f — OF  F POLL  TIME,  AND  EXHAUST  TAIL  IFNGTH. 

T A t-  s c values  will  default  as  foliows* 

V A 1 =•  PUNWAY-ARP  I VAL  INITIAL  VFLCXITY  = lAb  MILFS/HOUB. 

VA?  = RUNWAY -A»P I VAL  FINAL  VELOCITY  = 25  MILES/HCUP. 

V(  1 = FIINWAY-DEPAPTUPE  INIUAL  VFLOCITY=  0 MTLFS/HDUP. 

V L 2 - P UNWA  Y-Dfc  PAP  T UP  F FINAL  VELOCITY  = IfiO  MTLFS/HCUP. 

T I MF  = T A K F- OF  F MOIL  T I MF  = .01111  HOUPS. 

T/Il  “ FXAUST  TAIL  IFNGTH  = .8^23  MUES. 


ICI..PAUI  Tldl.Nf  .0) 

l F < Ai  I INPUT , 1 255)  VA  1 , VA? , VD1 , VL2, T IMF ,T AIL 
«R(  TF  If  111  PUT  ,9721  ) V A 1 , VA?,  VU1  ,Vr2.TIMF  .TAIL 
TFLTAIL.LE.CJ  T A 1 1 * 20./ 5280. 

VA  K * V A 1 * V A 1 
V V'  2 - V A?  « V A 2 


00 000164 
00000165 
00000167 
0030011  7 
00000168 
00000169 
00000170 
00000171 
00000172 
00000173 
000001 75 
00000175 
00030176 
00300177 
0000017F. 
GO  300 1 79 
OOOOOIPD 
OOOOOlfc 1 
0000018? 
00000163 
00000185 
00000165 
00000186 
00000187 
00000  1 88 
00000189 
00000190 
0000019  1 
00000102 
00000193 
00000195 
C000C195 
00300196 
00000197 
00000198 
00300199 
00000200 
00000201 
0030020? 
00000203 
00000205 
00  330206 
003002C6 
C0000207 
C00002C  8 

00000209 
00000210 
00000211 
00000212 
CO  300213 
00000215 
000002 l 6 

00300211 

00000217 

00000216 


i r 


XHIS  r’AGE  IS  BE-ST  QUALITY  PRACIICAW^ 
FMOM  Cur  Y i UiuXiSHM)  XO  DDfl  . -» 


r 


c-s 


c 

r 


VC  12  = VDI*VD1 
VC22  = VC2*VP2 


COMPUTE  RUNWAY  EMISSIONS 


>2 


DC  42  l * 1 » NPT 
J = IJ(U 
FMLMJ)  = 0.0 
fMTMJI  = 0.0 
TERM! J)  =0.0 
FMTxHJ  > = 0.0 

fmtxcmj  i = o.o 

«-  K ic  l J » = o.o 
FAPUtJl  * C.0 
DC  « K = l.NALT 
T NF  NO  = N AC  (K)  * NG  INIM 

fPLMJJ  = T7FNG*FINPGIK)*EM11INGMK  J,3,J)  •*■  FMLNIJI 
CXIMJI  = T NCNG*F  TKOF  ( K ) * EM  I I I NGN  I K ) , 4, J ) F FMTF(J) 


nr  B412  I At  C = 1,2 
PC  b2  1 « * l.NSR 
J S = 2RECPIMI/5280. 

CALL  l INF(XPFCPCM),VRFCP(M) ,ZK,X1 , Y 1 , Z l , X2 , Y2 ,2  2 ,C UP  » 
C SOI  FfRMATIlHO, I ItlX, 13(69.3, IX)  I 
rrupp  = LUP*r.(NVF/rtS 
nc  b21  LL  = 1,NPT 
J = ! J< l L ) 
cMSl  = FVLN(J) 

IM  lALC.tO.2)  c 3SL  = EMTF(J) 

PFCLAT(3,J,P»  = RFC.  0AT(3,J,M)  ♦ COUPP*FMSL 
#fCrMU,J,F)  = kECPAT(4,  J,M|  ♦ COUPP*FMSl 
c2 1 CCNTINUF 
rM2  CCNTINUF 


l TAXIWAYS  AS  LtNE  SOURCES. 

(..  I M PU!  AIRCRAFT  TAXI  S P F F C S l M I / HR » IN  THE  TERMINAL  AREA,  CN  THE 
(...  INPCUND  T AX  I WAY , ANO  C N T HR  OUTBOUND  TAXIWAY  RE  SPEC T I V 1 1 Y . 
IMDFAULT  ( 4 I .Nfc .0 1 
1»FAC( INPUT , 25 M FT AX  I (I) , 1 = 1 ,NALT) 

IF(UFAULT( 10) .NE.OI 
1 «FAG( INPUT ,25) ( FT  XI 1 (I) , I =1 ,NAl TI 
I F IFF  AUl  T I 1 1 ) .NF.O) 

IFF  At  < INPUT,  2SX  FIX  mm  , I = l,NACT) 

WRITE (CUT  PUT ,592) (F  TAXI (I) , 1 = 1 ,NACT) 

W°ITF( OUTPUT,  59 3)<FTXIl(  I >,!  = 1,NACT) 

« R l T F.  ( ( U T P U T , 59  « I (FTXICit  I 1*1=1, NACT) 

C 

C..  IJPLT  INITIAL  0IM7NSILINS  o TAXIWAY  WITTH  ANr  HFIGHT. 

( 

If(L'FAULT(I2I.Nr.C) 

IRFAGt INPUT, 3231  USTW.HTW 
,"5Th  * OSTW/2.4 
hT*  = HTW/1.2 
R l =0 


00  >00219 
C000022C 
00000221 
0010022? 
00000223 
00000224 
00000225 
00100226 
00000227 
0000022  P 
00000229 
00000230 
00100231 
C000023? 
00000233 
00000234 
00000235 
000002  3f 
00000237 
00000238 
00000239 
00000240 
00000241 
00000242 
00000243 
00000244 
00110245 
00030246 
00000247 
0000024  H 
00000249 
00000250 
00100251 
0000025? 
00010253 
00010254 
00000255 
0000025 6 
00000257 
00000258 
00000259 
000002CD 
000002C1 
0000026  ? 
0000026  3 
00000264 
00000265 
00000266 
00000267 
00000268 
00000269 
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3 I 


W R I TF  (CUT  PUT  , lM(FMir>(  IJ(J))«J*l,NPT) 

FCRMAT | IM*, TBO, •••»/» IX.  ••  TERM  APIA  FNG. 

WM  IF  (CUT  PUT  ,?*)!  (>  Ml  T(  l JU  > , l 1 , J«1  ,NPT  1 

FCRMATllH*, TBO, ••*,/, IX, ••  A 1 R(  RAFT  ARt  A • ,9X,  M F 10.  1) 
wBITFK  UT  PUT  ,301  ( t K IT! (J < J I , 2 I , J« l , MPT » 

FORMAT ( |M», TbO, *•*,/• IX, ••  NON-AIRCPAFT  AP F A • , 5X  , 5 |F1 0 . 1 ) 

FrwMAT(lll»,TRO,  •♦*,/,  IX,**  SEPVICE  V FFl  | Cl  F • , 7 X , 51  F 1 0 . 11  I 
IF  INAVF  .f  0.01  (.11  TU  1001 
* P I T t U. . 3 3 1 

FORMAT (1H*,T80, *•*,/, IX,**  ACCF9S  VFH1CIF*I 
r>r  m M*l,NAVR 

h»|TMl  UTPUT  ,36)  N,  ( A13F  I K,  1 J(  Jll  , J«1  ,NPTI 
FORMAT  1 1H»  , TBO,  ••*,  / , IX,  •••  , 16X,M,  4X.61F 10.1  I 1 

C r\T I\UF 

,«I1I  ((  UTPUT , 3d  I 

f CkM ATI  1M, , TbC,  ' * • ,/, IX ,791  • *•  1 ,//) 

WO  1 TF  |(  UTPUT  , 1000  1 
» ri’MAT  ( 1M1  I 
r CM  1 Not 
•T  U.  M1.F5P 

wp|  TF  11  UTPUT  , 81  1 NRUSIMM»XPFCP!N)»YRFCP|N1»(K‘AUTI!J(JII,J*1,NPT)  000006P  l 
FrwMATl  IX, 7V( • •• I,/,1X,,*,,T80,***,/,  000)066? 

1 1 X « • * • , 9 X , • Pf  LI JT  ANT  CIINC  FNTRAT  U N9  : F/CM  = M!CP0GFAM<  PIP  ( URIC  M F 0000  <T6H  3 
- T ep  • .THQ,*** ,/, 1X,***TH0,  COO 006 8 A 

2,**,/,lX,**'»l9X,*REClPTCH**  1 4,  TbO, •**,/,  IX,  •*•  ,T  80,  •**,/,  IX  00000686 

, 1VX, *f CCROINA TF s 1 * ,F7.3,*,* ,F7.?,* 1* ,T80,*** ,/,?(!*, •*• ,T800l  >00666 
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00000677 
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..*•*,/  1 , IX,  •* *,  31X,6IAb,?X)  I 
WO|TF (CUTOUT ,e3l  (P ter  ATI  1, 1 Jll  > , N)  ,l*l,KTT  | 

FC wMAT 1 IH, , TbO, ***./, IX ,* * AIRCRAFT  1 I N F *,11X,5(F IC.3) 
we  1 If  (CUT  PUT  ,84)  (PFCDATl  2,  IJIl  ),M,l*l  ,NPT) 

► TPMAT l IMF , TBO, *•*,/, IX, ••  AIRCRAFT  AFFA  • , 1 1 X , S 1 f 1 0 . 3 1 ) 

wPITF  ((  UTPUT,  82  MPFC  CAT  (3,  I Jll  1 ,N  I , l * 1 , NPT  I 

1 CKF'AT  1 1M6.T80, ••*,/, IX,**  N(N-AIRCRAFT  APfA-ltNFS',2X,5(Fl0.3ll 
wPI  Tt  (I  UTPUT  ,861  (PFCOAT  (4,  IJ(l  I ,N1  .1  *1  , NF'T  1 
f CRMAT1 IH6, T80, **',/, IX, •*,,12X,* TOTAL* ,9X,6(F10. 311 
wR I Tf (OUTPUT , 38) 


LCfT  I NUr 
. . INPUT 


F ( F MA T S . • . 
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F C PM A T 1 l2,?Fh.u» 

FfKMAT ( 16X, H 61 
FORMAT  1 SX,  IMS) 

' rPMAT  UX.231  3) 

FORMAT < 6 X , I 6 * 6 • 0 I 
» rPMATI16X,*F 10.0) 

I CRM  A T I l 9 X , ?F 10 . 0 ) 

F CpmAT 1 16X,6F10.0) 
frFMAT(10X,(F10.0,F8.C,l2l 
F f PM A T 1 IX, A6.5X, 10F5.1  I 
F ( PMAT1 16X.4F IC.O) 
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f CASrilNE  ENGINE  POWEPEr  MOTOP  VEHICLE  EMISSIONS  ARE  CCMPUTFr  ON 
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CAN  “V  (JStr  ME  SMALl  -ANf.lf  GASF  hMFPr  LINE  fS  FL/PTHFP  SEGMENTED. 
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Xl*X 

Y l*v 
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Y 2*YB 
/ 2 * Z B 

X«0.5*< Xl»X2) 

Y»n.s*(Y i ♦ y r » 
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i i X * XX 

Y * YA 
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tf  < ISAC  .FC.  1 t nwoi*  xp  -xi 
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! f I X 1 . f <-.*2  ) 01  TO  44 

Z=Zl*lA-Xl)*(Z2-/ll/IX2-Xl) 
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S I GV  = S I GZ  1 J c T At  ,Th4V) 
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l*  A 1 = - * A»f,l  »APG21  » / I 1.4l42*SIGH*SIGV*RArG> 

HA*  = - ( ARG1*  ARG2  2 ) /(  1 .-VI  42*S)GEl*SIGV*P  AEG) 

APUY  l = Al *hA  l 
APGY2=AL*hA2 

I 1=  t YE-Y  l )**.  /LENH**2*  ( Z*-  - Z 1 ) **2/0t  NZ**2 
C 2»(  YE-Y1  )**2/DENEK  *2E(ZRFZl)**2/DfNZ**2 
AR&Z  > = EAl**2-Cl 
ARGZ2r8A2**2-C2 

44 t l F ( At-GZl  .11  .-10.  . ANP.APGZ2.LT  .-10.  ) GO  TP  500 
I f ( ARGZl .1 T.-IO. )G0T0  2411 
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V A --  SON  T ( F re  T A | ♦ VS 

V * SCPTlPCPTei  ♦ VS 

YA  = SOP  T ( WS2  * VA  * I V A ♦ WSCII 
YP  * SO0  T ( N S2  ♦ VI'  * I VH  ♦ NSC  1 ) 

\o.,  = (Yh  ♦ VO  *■  NSr/2.»  / ( YA  ♦ VA  ♦ WSC/2.1 
;i  * YW  - YA  - NSC / 2 • * A LOG I AFG) 

•:i  = K / TAIL  * (FXSI  / SP  CL  / A) 
f'FTUPN 
= NP 


00  TO  l A 1 3 
0000  l A ( A 
0000141S 
00001416 
00001417 
‘0000 1 A 1 6 
00001419 
-00001420 
00001421 
00001422 
00)01423 
00001424 
0000 1425 
00001426 
00001427 
00001428 
00001429 
00001430 
00001431 
00001432 
00001433 
0 0001434 
0000 14’S 
00001436 
00001437 
00001438 
00001439 
000014*0 
00001441 
00001442 
00001443 
00001444 
0000 144  S 
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F GnC  T JC  N Slf.GKt  TM( |*JR  ) OMCIUi 

CIHMON/STA6  II  /NiTb,  Sll.XC),  SIG2C  OOOOl  A A 7 

9lGGX*SlGYlKST8,THUUP  l-SIGXP  00001446 

> f ruKf  oooo 1444 

CND  00001490 


M. Ni"  t 1 1 N S I GG/(  I MOUH  I 00  7014*1 

l i 88PN/9TA8U  /NS TP.  SIGXO.SIGZP  0Cr0l49? 

9 1<.G2*S  !G2 I X9T8, TMOUB 1-SIG2P  00301493 

► r T U6N  0 000  14  54 

197  0000  l 4 5 *• 


* GNC  Til.  S1GYI  J . TMIIUP  ) 

C.OMCIN/M  liN/v.'-A Vfc 
Ml  NS  It  N A < 71  , *<  71 
•'TMFNSli  A f I'.I.JIM 

rAIA  l /I  99.  . ICO.  ,68.  ,90.  ,34./ , EV. 91  ,.««*,  .43,  .40,  .93/ 

” AT  A A/2.191  l ,1  .949-*,  1 ,Ut06,  .6  846  9,  .99361 , . 69  36  t,.993fco/ 

C A T A ts/. 6732  6, . d82ol. .69031.. 86  866, ,8<  138,.  pc  138,.  89  n6/ 
r / r A < (.  N V / • OOOt  213  71/ 

t t t * r ♦ • ^«««r#*«%«*»***««#****«*»*******»*f.  •***•♦<**#♦ 

rr^^yTt  HOPIZrNTAl  riSPt»>SlON  frFFMCirNT. 

i t « *e  *«*«*t***t  •*#*>*»«*«•**<***«*•*•**  **•«»•%***«-?***•******< 

( J*l*2»3»4»9  APf  ClASSF*  8,C,l,r,F 

( X«  4f.f  SIGXV  APl  IN  MIlFS 

T9f(  -THUURo’tOO. 

i 

x I •«  Y*  l A 1 J | * 1 SEC  * • 8 ( J)  I ” C.I  N V * • H 7 
f ».S.\VC  19  Mf  AN  V.  I NO  SOCK  IIP  TO  Tl«f  T. 

>X  = W SA Vt  * T Of  UK* 1 .609 
I J 5 9 A XO I J-l  , 1 ) 

SlGXY*C 1JJ|»1XX**DIJJ) >*.000889 
l .C00tl9« J 1/ I6f9. »*(t**.2» 

9 ll.Y  ,/MAXl  ( SIGY  , S IGXY  I 
f 

(.  PISPlPilllN  UMFIC.1FMTS  I AS  E 0 ON  TllRM  E W(.PKF*n(lK. 

f bl-FM  Till  9 1 f.»A  s GIVEN  IN  THI  P.OPKPCCK  P|R 

' Tf  ti  8 1 MIT  F S.A8PI  INGS  APT  SCAIEO  TO  1 KUP 

F A V * P A F E S «Y  A F A T.  IP  or  1.43 

l‘’Ui<l, 

r a > 


0000  1496 
0C7014*  7 
00001466 
00  101499 
00  70  |v  C 
00)0141  1 
CO00 146? 
roioi4<  3 
c nor  ip. a 
0000141' 
0070  14( 6 
000014/  7 
000014/ 6 
00001469 
C03014  70 
0000  14?  1 
0070147? 
0030 14  7 3 
00001474 
00^0147* 

0000  14  7/ 

0000 1 A 7 3 
00001- 
00  JO’.  A 
00 70  1- 
O00i> l a : 
00  70  1 - 1 * 
00701463 

0000 1464 
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F UNC  I IC  N Sir./IJ.THfHJR) 

C C*M(  N/*njN/WSAVF 

( I I SPEf-  1 S I f N CPM-.  FASFt'  ON  TUPNFPS  hCPKBnt* 

r j=i ,2, 3. a, 5 ape  classfs  B,r,n,p,p 

l X ANP  SIC./  ABF  IN  MILTS 

f'lNlNSU  N C 13,51 ,01 3.51 
•UAL  TIME  ( 71  , A ( 7,6)  ,B(  7, A) 

C 

TATA  1 I ML /Q. , 300. , 1000. , 3C 00 . , 1 OOOO ., 30000 . , 17/000./ 

PAT  A A/.  I 71  //,./  760  8,  .4  l /IS,  .519/1  , . 5096  3,  .A  7/3  9,.  5/1  AO, 

1 . 1 101  2,.  3955 3,. 4 1/19,. 57145,. 76485,.  7 1 936 , . B8 88/ , 

? .0133b,. 16640, .A  12 19, l.CBl  3,  1. 946 7,/. 39 3 t, 1.6(77, 

I . C 1336, .1 6640, .41/19, 2. 2830,/. S« 50, 3. 8684,6.745/. 

3 .013‘'8,.166AC,  .A  1/19,/.  3 33  3,5.7990, 16.(97, /O.  673  , 

3 .013  38,.  1664C,  .41/19,5.6  801, 14.599,64.577,54. 149/ 

AT  A.  8 / 1./058, 1.057/,. 9/365,.  84130,.  796  8 9,  . 76  308  ,.  69  839  , 

1 1 .ZB  64,. 99/ 7 5, .92365, .b/449, .7/571, .6508 /,. 60486 , 

/ 1.59//.1 . I 195, .5/36  5,.  73/1  7, .5904  7, .51  700. .49  583, 

1 1.55?/, 1.1195, .9/ 365,. 63 88 3,. *3706 ,.4  5686,. 336 77, 

•.  1 . 89//,  1 . 1 195,  .9/36  5,  .6  3646,  .4649  7,  . Z962  1 , .21517, 

8 1 .89/?,  1 . 1 1 95, .9/36  5, .55016, .37541,.  16  60 7, .1/177/ 

f A 1 A C/1  10., 110. , 11C. , 60. , 60. , 60. , 33. ,33., 40. ,/l. 5. /I. 5, 36. ,14., 
1 14. ,23. 5/ 

AU  P/1. ,1.09, l.o*,. 92, .9/, .92, .80, .61, .53, .70, .56, .35, .78, .53, 

1 . 30/ 

rATA  rp\V/.OPOo/1371/ 

( CPMPuU  VFCTICAl  P!SPrrS!nN  COEFFICIENT. 

T S£C  *THflUP*  3600. 

CC  10  N = 2 , 7 

I F( TSFC ,Lf . T fMf (N)IGO  TP  20 
10  C t NT INUF 

ILlN.r,  T.7I  N = 7 

C T I Mf  OF  TF  AVFL  S^OULO  Ht  LF5S  THAN  172. F03  SFC  APPPPX.  2 OAYS 
>r>  f t NT  1 NUF 
N = N-1 

S I .AM  A ( J,M*TS-fC**b(  J.NI  l*CC>NV 
C T.)F  N F 8 TFST  FoUA  T I IN 

( 5 11,2*  I 1.17518*1  S=C**.7/Cl/)*CrNV 

XX  = hSAV-«-THri!S*  1.609 

ru 

1 F ( XX. (T .1.1  1=/ 

IMXX.OT.IO.)  I - 3 
J J^AXOI  J-l  , 1 > 

4 T - , T / = (f ( l , JJ  |»x<*»1|| , JJ) >*.000714 

( ,r  I0714=( 1/I6r9>*(20/10)**.2,  ASSUMING  THAT  THE  VFSTICAL  PISPEP- 

( s 1 1 A<  Cc  f 6 P I f I FN  T IS  INSEA  SI  T|  VF  TP  SAMPLING  TIMFS  OVTF  20  MTN. 

1 ,/^AMAXl<riG2,SIGT/l 

- r T IJK  *. 

•:  An 


00001485 
00001486 
00001487 
00001488 
00001485 
00001490 
0000 1491 
0000149/ 
00001493 
00001494 
00001495 
00001496 
00001497 
00001498 
00001499 
00001500 
00001501 
0000150/ 
00001503 
00001504 
00001505 
00001506 
0030 150 T 
0000 l 5^ 8 
00001509 
00001510 
*00001511 
0070151/ 
*00001513 
00001514 
00001515 
00001516 
00001517 
00001518 
00001515 
00001520 
00001521 
000015/2 
000015/3 
00001524 
000015/*- 
0000 15/ f 
000015/7 
000015/8 
000015/5 
00001530 
C0001531 
0000153/ 
00001533 
00001534 
00001535 


L 


C-35 


FUNCTION  TPAF  ( NS  ,KS  T AF  , HE  F A , NP  Ft  AG  , I RAC  K I 
f CKMCN/OF  IT  A/CP  LY.OFLZ 

( tMMPK/ior /rw,c * 

(TMMON/MPT/WS.WD,  JST  AB.FIMXl  ,TFMP,CL,SUCC)Y,  SUOOZ 
CCMMON/PL  3/ X S( 5 I » YS ( 5 1 ,STKH(5)  * W l T 1 ( 5 ) 
CCMMuN/PFCPT /HPFCPA.HTAEPO, ZKECPC 
CCMMCN/SPUEC /IX, TY, T2 
CCMMnN/XTPAN/WSMD.NCALM, S02PI 
CCMMUA/nDUN/hSAVE 


0000  I 536 
00001537 
0000 1536 
000015351 
00001540 
00001541 
00001542 
00001543 
00001544 


f 

( 

f 


COMPUTE  COUPLING  COEFFICIENT  AT  PECEFTOP  CUE  TO  POINT  ANE  AR F A 00001546 
SOUFCE.  00001547 


r IFlKSTAe.GT.OJHAVF  PLUMF  INITIALLY  ABOVE  LID. 

2S=STKHINS)/52B0. 

HfEF=UFFA/52BC. 

» I T = W I T 1 1 NS  ) 

JSTA8T=JSTA8 
wSAVE  = 8 SMC 
HRtCPT=HRf CPA/5280. 

2KEFP  = F“<FCPT 
A001=0. 

AC  02- 0 . 

ACi)a*0. 

ACD4=C . 

Ai'Os-O. 

A0D6=C. 

*>  J wO  - 0 • 

XL  = CL 

C 

I FIKSTAB.GT ,0»G0  TO  121 
taio=HPixi/E2ec. 
r.n  TO  140 
121  JSTAEsT  = 4 

IFIKSTA8.E(j.2>  GOTO  1 30 
JSTABT=5 
130  CONTINUE 

HL 10=1 000 . / 5280 . 

X L = 30 . / 52  HO. 

140  CC  NT  I NUF 

( IF  PLUME  p,  ELCw  LIG  ANI  PECEPTCR  AbCVF,  TPAN=C. 

I F (UP  ECPT.GT  .E-L  I0.AN0.KSTAB.E0.0IG0TC76 
I F ( * S A VF  .Lf.  WJWGI  GO  TO  9596 


0000 1549 
00001550 
00001551 
00001552 
00001553 
00001554 
00001555 
00001556 
00001557 
00001558 
00001555 
00001560 
00001561 
00001562 
00001563 
COOC 1564 
00001565 
00001566 
00001567 
00001568 
00001569 
00001570 
00001571 
00001572 
00001573 
00001574 
0000 1 5 7 c 
0000  l •>  76 
00001577 
00CC157,- 


PMN  = APSICW) 


00001579 


T M I N = Uw/  i»S  A VF 
706  C TNT  I NUF 

T T *T  M I N»T Y 
IF  ITT .1  t.C. 1 GOTO  76 
IF (NPFL AG.NF.01  GO  TO  143 
T*1N  = 0. 

IF  (l)W  .OF.  WIT!  TM I N = LOW  - WIT)  / WEAVE 
I F ( IF.f  CK.FO.O)  GO  TC  131 
IFIEWFWl  I -Of  L Y)  132,132,123 


00001580 
0030 1581 
OCOO 1582 
0000 1583 
0000 I5f 4 
0000  158  5 
00001586 
00001587 
00001588 
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I!?  TMIN*C. 

Of  TO  131 

13  3 TM|N*(UW»Wl  T-rUYl/WSAVF 
131  TT  * TM)N  ♦ TV 

IT  (TT.lt  ..01  ITT*. 01 
<Y*SILY ( JSTAHT.TT 1 
XT  P * Vi  T T ♦ r* 

If  ( XT  P .0*.  2.  *H  1 T I C.(  TT  IA2 
Pk  * XPP  / 2. 

1A2  T IT  * TZ  ♦ f W / WSAVF 
s / » S I OZ  < jSTaht, r r t I 
‘ P*DMIM/SV 
V XPl «-.S*f  P*f  c 
I F ( t XPP  . I T , - 2C.  1 Gl'TP  76 
f CO*EXP( FXPT 1 
l F | Oh  .GT  • 2 • * XI  KIOTO  153 
riM.M-P.2H  H H53*WSA Vf  *SY*5Z 
7 1 * HFfP-ZKFFP 
E XP/l*(-.s« (ZI/SZ)*(Z1ZSZ>» 

If  I P X P Z 1 .1  T .-20.  lOliTO  76 
• l*k  XP(I XPZ  1 I 

AT 01 *F 1 *f  or/r  FNOM 
IF  (Of  FF  .GT.r.Ol  GO  TO  171 
xru2*Arm 
GCM  172 
1 71  CfM  INUt 

1 7 * HEFF  ♦ ZKC  CP 
f XPZ2*I-.5*I22/S2)*IZ2/S/  I) 

I F (f XPZ2.1  T .-20. l&OTC  M 
F2*f XP(f XPZ21 
tC"'2*APri*F?/Fl 

172  rrNTiMjr 

I FlhfFF.GT.Hl  I0/2..0F  .SZ.C.F.(Hlir-FEFFl/2.2l  GO  TO  If 
GC  TO  tl 

10  UNTINUf 

IFlhfFF.LT.MUri  GO  TO  1 7<* 

ATG3*  ACM 
AruA*ACC2 
GC  T ( 1 7 3 

17a  CfMIMJF 

Z?*  ( 2.*Ml It-HFFF l-ZKcFp 

I XV»/3=(-.«>*(Z3/i/J*1^3/S21  1 
I f (FXPZ3.1T.-20. (GOTO  t 1 
■ ? *£  XP  ( I XPZ3I 

7ri)A*..l01*F3/fi 
/A*  Z3*2.*ZXFFp 

I %\>i  A=-.5«(  ZA/S/I*(  ZA/r,Z  1 
|F11XPZ'«.IT  .-20.  )G(1T0  61 
XPZA  1 

Af  OA*  A C IO l *F  A/F  1 

If  (Hf f F .GT.O.OI  CP  TO  173 

t rGWl.'t>3 

A C r><  = A C I A 
G r II,  u 

173  f.  CM  | v til 


0000 1 5P9 
00  30 1 590 
0000  1 591 
00001592 
0000  13113 
0000 l 5 VA 
00001595 
00001596 
00001597 
00001598 
00001599 
00001600 
00001601 
00001602 
00001603 
0000 160A 
00001605 
00001606 
00001607 
00001608 
00001609 
00001610 
00001611 
00001612 
00001613 
00  30 16 1 A 
00001615 
000016 16 
00001617 
00001618 
0030161* 
00001620 
00001621 
00001622 
00001623 
0000 162 A 
00001625 
00001626 
0000162T 
00001628 
00001629 
00001630 
00301631 
00001632 
00001633 
00001 63 A 
00001635 
00001636 
00001637 
00001638 
000016.39 
000016A0 
000016A1 
000016A2 
003016A3 


*»IS 


corf Is  S£sr  aUAr 
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Z5-Z242.MHI.  IC-ZKFFP) 

00001644 

r xp/5*-.->*iz5/sz)*<  zo/sz» 

00001645 

IMF  xfz5.it  . -20.  ir.nrn  m 

00)01646 

F *«t  XPlf  XPZ5) 

00001647 

AC05*ACPl*F  5/FI 

00)01648 

Zf «Z5*2.*ZKFFP 

00001644 

FXPZ6=-.5*1Z6/SZ)4( Z6/F  7 ) 

00)01650 

lF(FXFZ6.LT.-20.)GQTn  61 

00001651 

f e*t  XHF  XPZt  1 

00001652 

Aro6*AtCl*l t/F 1 

00001653 

r.c  Tf  el 

00)0  1654 

1S3 

PFM3M*  SC2P  1 •w5AVF*SY 

00001655 

ACul*tDP/<rFAOM«Hl IP) 

COOO  1656 

61 

rr\T  imjf 

00001657 

THAN- AT  PI ♦ATT 2* A0D3*Arr4*A0n5  + A0C6 

00001658 

I-  F TUPN 

00001655 

143 

1 F I InAlK.tO. 1 ) 00  rn  144 

0000 1660 

3 2*(Cki4«tIT)/  VSAVE*-TZ 

00001661 

1F1FW-WIT)  1 45* 145,1 46 

00)01662 

14‘ 

T 1 = t Z 

00001663 

T T = T Y 

00001664 

GC  TC:  145 

0000166  5 

l 46 

T IMl  Yw-m  1 T ) ✓ WSAVFfTZ 

00001666. 

r T*(n»-h  m / *»s a Yf ♦ t y 

00001667 

( r TC  144 

00  )0166e 

144 

r?=Xl/  WSAVFtTZ 

00)0166  5 

I F ( otutw  1 T-PF  1 Y)  147.147.14H 

00001670 

147 

T 1 *1  Z 

00001671 

n = TY 

00001672 

3C  TC  144 

00001673 

1 4t 

T 1 = 1CH*W  1T-THLY  1 / oSAV^TZ 

00001674 

r T-tuk*-y  1 T-r  t LY  ) / d S A VF  * T Y 

COOO  1675 

14! 

IF(Tl.FQ.lZ)  GC  in  iso 

00)01676 

sz i = sigz i js  taht , r n 

00001677 

cr  Tt.  151 

00)01678 

15L 

szi=  mz/2.4 

00001675 

151 

SZ2  = F 1GZ ( JSTAPT.T2) 

00)01680 

IF(  TT.Ff3.TY)  Gf>  T 0 152 

00C01681 

5y=sigyi  jstapt.tt) 

00001682 

GC  TO  155 

00001683 

15? 

S Y = f-  £ L Y/2.4 

00001684 

1 "»6 

A«fiY=-<C  w/SY  )*<  2/2. 

000016-  •> 

t F ( APOY.l  T.-2C.  1 (.  ) Tf.l  74 

0000  16(  . 

h=A1CG(SZ1/SZ2)/ALHC(T1/T2| 

OOOOlM  7 

1 54 

r X Z l = C . 

00)016-  , 

FXZ?=0. 

00)0l6t-' 

ACZl=-<  (ZS-HFFC.PT  I /SZ2  l*o?/2. 

00C016SC 

A CZ  2 = - l ( ZS  + F-F  FC  PT  »/SZ?l  **2/2. 

00)01641 

IF  (At.tl.LT.-20)  GC  TC  15  7 

00)0169? 

) XZ 1 = 1 XP<  ACZ  1 ) 

00)0164? 

Iff 

1 F ( Af.ZZ.l  T . — 201  GC  TC.  158 

00)01654 

*■  xz*  = t XP  1 A(>Z  2 ) 

00001645 

15H 

f CM  IMJF 

00)0164< 

|F(|EXZ1»FXZ2).LF.0.)  GO  TO  76 

00001657 

F X Y = 1 XP ( APGY ) 

0000164« 
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ihau;>im-i.  ).u  ..on  or  to  2 

OOOOtfcSS 

M Mk>IT2/S??>Tl/Sn)/li.-t*l 

00001700 

Cl  Til  l 

00001701 

> 

l rNX»Tl*All‘C(T2/Tl»/SZl 

00101702 

tpak»h»v*u  x/i*U22»*fCNx/u  sc2Pt  **2»*sv*dfi  vi 

0000170) 

l i*Al  *»C 

00001704 

wf  TUFA 

0000170) 

lb 

T»AA*0. 

00001704 

tPAl  A*(J 

00001707 

uf  TUHN 

0000 1708 

Cum  i mjf 

0000 l 70S 

l rt A(  K--0 

00001710 

uMUfK 

00001711 

pAP 

00001712 

